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FOREWORD 


Structure,  propulsion,  and  guidance  of  new  or  improved  weapons 
delivery  systems  are  dependent  in  crucial  areas  of  design  on  the 
availability  of  accurate  thermodynamic  data.     Data  on  high-temperature 
materials,  new  rocket  propellant  ingredients,  and  combustion  products 
(including  exhaust  ions)  are,   in  many  cases,   lacking  or  unreliable. 
A  broad  integrated  research  program  at  the  National  Bureau  of  Standards 
has  supplied  new  or  more  reliable  thermodynamic  properties  essential 
in  several  major  phases  of  current  propulsion  development  and  application. 
Measured  were  compounds  of  those  several  chemical  elements  important 
in  efficient  propulsion  fuels;  those  substances  most  affecting  ion 
concentrations  in  such  advanced  propulsion  concepts  as  ion  propulsion; 
and  the  transition  and  other  refractory  metals  (and  their  pertinent 
compounds)  which  may  be  suitable  as  construction  materials  for  rocket 
motors,   rocket  nozzles,  and  nose  cones  that  will  be  durable  under 
extreme  conditions  of  high  temperature  and  corrosive  environment.  The 
properties  determined  extend  in  temperature  up  to  6000  degrees 
Kelvin.     The  principal  research  activities  were  experimental,  and 
involved  developing  new  measurement  techniques  and  apparatus,  as  well 
as  measuring  heats  of  reaction,  of  fusion,  and  of  vaporization; 
specific  heats;  equilibria  involving  gases;  several  properties  from  fast 
processes  at  very  high  temperatures;  spectra  of  the  infrared,  matrix- 
isolation,  microwave,  and  electronic  types;  and  mass  spectra.  Some 
of  these  techniques,  by  relating  thermodynamic  properties  to  molecular 
or  crystal  structures,  make  it  possible  to  tabulate  reliably  these 
properties  over  far  wider  ranges  of  temperature  and  pressure  than  those 
actually  employed  in  the  basic  investigations.     Additional  research 
activities  of  the  program  have  involved  the  critical  review  of  published 
chemical-thermodynamic  (and  some  chemical-kinetic)  data,  and  the 
generation  of  new  thermo chemical  tables  important  in  current  chemical- 
laser  research. 
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ABSTRACT 


This  report  (covering  July  1973 -June  1974)  gives  accounts  on 
the  determination  of  the  thermo physical  properties  of  conductors 
by  high  speed  techniques,  studies  of  reaction  intermediates  of 
possible  interest  to  chemical  laser  systems  by  the  methods  of  matrix 
isolation  infrared  spectroscopy,   preliminary  report  on  the  reaction 
kinetics  of  Ba(g)  with  vibrational ly  excited  ozone,  a  survey  of  the 
dissociation  energies  of  the  scandium  group  and  rare-earth  oxides, 
and  a  bibliography  on  spectroscopy  of  fluorides  and  oxides  belonging 
to  the  lanthanide  series. 

Using  a  subsecond -duration  transient  technique  the  specific 
heat,  electrical  resistivity,  and  hemispherical  total  emittance  were 
simultaneously  measured  over  the  temperature  range  1500-3000  K  for 
some  grades  of  graphite.     Similar  measurements  were  made  on  vanadium, 
and  zirconium  in  the  temperature  range  of  1500-2100  K.  Melting 
points  and  radiance  temperatures  (at  650  nm)  are  reported  for  zirconium 
and  molybdenum.     The  temperature  of  the  transition  from  the  o/-(3  phase 
of  zirconium  and  the  energy  difference  of  these  phases  has  also  been 
measured  using  the  subsecond  duration  transient  technique. 

The  products  of  the  reaction  of  Fe(g)  and  O2  have  been  investigated 
and  identified  using  the  methods  of  infrared  matrix  isolation 
spectroscopy.    A  preliminary  report  on  the  study  of  the  reaction  of 
Ba(g)  with  vibrationally  excited  ozone  is  also  presented. 

The  dissociation  energies  of  the  scandium  group  and  rare-earth 
gaseous  monoxides  are  evaluated  by  reviewing  the  literature  available  on 
Knudsen  effusion  rates  and  mass-spectrometric  data.    Criteria  are  discussed 
for  choosing  data  for  evaluation  of  dissociation  energies.    A  bibliography 
of  the  available  literature  on  the  spectroscopy  of  fluorides,  oxides, 
and  oxy fluorides  belonging  to  the  lanthanide  series  is  given. 


Charles  W.  Beckett 
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Chapter  1 

MEASUREMENTS  OF  HEAT  CAPACITY,  ELECTRICAL 
RESISTIVITY  AND  HEMISPHERICAL  TOTAL  EMITTANCE 
OF  TWO  GRADES  OF  GRAPHITE  IN  THE 
RANGE  1500  TO  3000  K  BY  A  PULSE 
HEATING  TECHNIQUE* 

A.  Cezairliyan  and  F.  Righini 

National  Bureau  of  Standards 
Washington,  D.  C.  20234,  U.S.A. 

Abstract 

Measurements  of  heat  capacity,  electrical  resistivity  and 
hemispherical  total  emittance  of  Poco  and  pyrolytic  graphites  in  the 
temperature  range  1500  to  3000  K  by  a  subsecond  duration  pulse  heating 
technique  are  described.     For  a  given  graphite  grade,  heat  capacities 
of  different  specimens  were  in  agreement  within  0.5%.     The  difference 
between  the  results  for  the  two  different  grades  was  about  1.8%;  the 
results  of  Poco  being  higher  than  those  of  pyrolytic.  Electrical 
resistivity  of  the  Poco  graphite  was  about  four  times  greater  than  that 
of  pyrolytic  graphite  (parallel  to  basal  planes).    Hemispherical  total 
emittance  of  Poco  graphite  was  almost  twice  that  of  pyrolytic  graphite. 
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1.  Introduction 

Disagreements  exist  between  the  limited  number  of  experimental 
results  for  the  properties  of  graphite  above  1500  K  reported  in  the 
literature.    Most  of  the  reported  measurements  were  performed  on 
different  grades  of  graphite,  which  complicates  their  evaluation  and 
poses  the  question  whether  the  differences  were  due  to  measurement 
errors  or  were  really  indicative  of  differences  in  the  graphite  grades. 

As  an  attempt  to  elucidate  this,  a  program  was  initiated  for  the 
systematic  and  accurate  measurement  of  selected  properties  of  various 
grades  of  graphite  at  high  temperatures.    The  results  of  the  first 
phase,  namely  the  measurement  of  heat  capacity  of  a  grade  of  graphite 
(Poco,  AXM-5Q)  ■■   in  the  temperature  range  1500  to  3000  K,  were  presented 
in  an  earlier  publication  [1].    The  objective  of  this  paper  is  to  report 
the  results  of  the  measurements  of  heat  capacity,  electrical  resistivity 
and  hemispherical  total  emittance  of  two  grades  of  graphite  (Poco  and 

pyrolytic)  in  the  range  1500  to  3000  K.    The  Poco  graphite  used  in  this 

-fro  too 

work  (Poco,  DFP-2)  is  differentiae   one.   used  in  the  earlier  study.  For 

a. 

completeness  and  for  ease  of  comparison,   the  results  of  the  earlier  study 
on  Poco  graphite  are  also  included  in  this  paper,. 

^Commercial  identification  of  graphites  is  made  in  this  paper  in  order 
to  specify  adequately  the  grade  of  the  graphite  specimens  whose  properties 
were  determined  experimentally.     In  no  case  does  such  identification 
imply  recommendation  or  endorsement  by  the  National  Bureau  of  Standards. 


2.  Method 

The  method  is  based  on  rapid  resistive  self-heating  of  the  tubular 
specimen  from  room  temperature  to  high  temperatures  (above  1500  K)  in 
less  t;han  one  second  by  the  passage  of  an  electrical  current  pulse 
through  it;  and  on  measuring,  with  millisecond  resolution,  such 
experimental  quantities  as  current  through  the  specimen,  potential  drop 
across  the  specimen,  and  specimen  temperature. 

The  current  through  the  specimen  is  determined  from  the  measure- 
ment of  the  potential  difference  across  a  standard  resistance  placed  in 
series  with  the  specimen.     The  potential  difference  across  the  middle 
one  third  of  the  specimen  is  measured  between  spring-loaded,  knife-edge 
probes.    The  specimen  temperature  is  measured  at  the  rate  of  1200  times 
per  second  with  a  high-speed  photoelectric  pyrometer  [2],    The  small 
hole  in  the  wall  at  the  middle  of  the  tubular  specimen  provides  an 
approximation  to  blackbody  conditions  for  optical  temperature  measure- 
ments.   All  these  data  are  recorded  with  a  digital  data  acquisition 
system,  which  has  a  time  resolution  of  0.4  millisecond  and  a  full-scale 
signal  resolution  of  one  part  in  8000. 

Details  regarding  the  construction  and  operation  of  the  measurement 
system,   the  methods  of  measuring  experimental  quantities,  and  other 
pertinent  information,  such  as  the  formulation  of  relations  for 
properties,  etc.  are  given  in  earlier  publications  [3,  4], 
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3.  Measurements 

The  new  measurements  reported  in  this  paper  were  performed  on 
one  Poco  graphite  (DFP-2)  specimen  and  two  pyro lytic  graphite  specimens. 
The  characteristics  of  these  specimens,  as  well  as  those  of  the  Poco 
graphite  (AXM-5Q)  specimen  used  in  the  earlier  study,  are  given  in 
Table  1.    In  the  rest  of  this  paper,  the  specimens  are  referred  to  by 
their  designated  numbers  as  listed  in  Table  1.    The  pyro lytic  graphite 
specimens  were  fabricated  in  a  way  to  allow  the  pulse  current  to  flow 
in  the  direction  parallel  to  the  basal  planes.    Before  the  start  of 
the  measurements,   the  specimens  were  heat  treated  by  subjecting  them 
to  15-20  heating  pulses  (up  to  2800  K) .    All  the  experiments  were 
conducted  with  the  specimen  in  a  vacuum  environment  at  approximately 
10"5  torr. 

To  study  possible  effects  that  might  be  attributable  to  the  rate 
with  which  the  specimen  heats,   two  ranges  of  heating  rates,  designated  as 
"fast"  and  "slow",  were  used  in  the  experiments  on  specimens-l,  2  and  3. 
To  optimize  the  operation  of  the  high-speed  pyrometer,  the  temperature 
interval  (1500  to  3000  K)  was  divided  into  six  ranges.    A  "fast"  and 
a  "slow"  experiment  were  performed  in  each  temperature  range.  These 
yielded  a  total  of  12  experiments  for  each  of  the  specimens-l,  2  and  3. 
In  the  case  of  specimen-4,  only  a  "slow"  heating  rate  was  used  yielding 
a  total  of  6  experiments  for  the  specimen. 
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Table  1 


Characteristics  of  the  Graphite  Specimens 


Specimen 

ru  l x  l y 

UC  1 1  o  L  L y 

Nominal  Dimensions,  mm 

Designation 

Grade 

-3 

g  •  cm 

length 

O.D. 

thickness 

Specimen-1 

Poco  (AXM-5Q)* 

99. 9+ 

1.755 

76 

6.3 

0.5 

Specimen-2 

Poco  (DFP-2) 

99.99+ 

1.791 

76 

6.3 

0.5 

Specimen-3 

*** 

Pyrolytic 

99.99+ 

2.19 

76 

7.2 

0.5 

Specimen -4 

•kick 

Pyrolytic 

99 . 99+ 

2.19 

76 

7.2 

0.5 

The  graphite  block  was  kindly  furnished  by  the  U.S.  Air  Force  Materials 
Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio. 

"k 

The  graphite  block  was  purchased  from  the  Poco  Graphite,  Inc. 

The  graphite  in  the  form  of  tubes  was  kindly  furnished  by  the  Union 
Carbide  Corporation  through  Dr.  A.  W.  Moore.     Deposition  temperature: 
2400  K. 
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To  optimize  the  operation  of  the  measurement  system,   the  heating 
rate  of  the  specimen  (in  both  "fast"  and  "slow"  experiments)  was  varied 
depending  on  the  desired  temperature  range  by  adjusting  the  value  of 
the  resistance  in  series  with  the  specimen.     Duration  of  the  current 
pulses  ranged  from  210  to  330  ms  for  the  "fast"  experiments,  and  280  to 
520  ms  for  the  "slow  experiments.    Other  typical  operational 
characteristics  of  the  measurement  system  during  the  experiments  on 
the  two  different-grade  graphite  specimens,   such  as  specimen  heating 
rate,  and  heat  loss  from  specimen  by  thermal  radiation,  are  summarized 
in  Table  2. 

Optical  checks  performed  on  the  experiment  chamber  window  after 
the  pulse  experiments  did  not  indicate  any  graphite  deposition.  Also, 
weight  measurements  before  and  after  the  entire  set  of  pulse  experiments 
did  not  show  any  change  in  specimen  weight. 

4.     Experimental  Results 

This  section  presents  the  properties  determined  from  the  measured 
quantities.    All  values  are  based  on  the  International  Practical 
Temperature  Scale  of  1968  [5].    In  all  computations,   the  geometrical 
quantities  are  based  on  their  room  temperature  (298  K)  dimensions. 

4.1.    Heat  Capacity 

The  heat  capacity  data  were  determined  from  the  measurements  of 
current,  voltage,  and  temperature  obtained  during  the  heating  period. 
A  correction  for  heat  loss  due  to  thermal  radiation  was  made  using  the 
results  on  hemispherical  total  emittance. 
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Table  2 

Operational  Characteristics  of  the  Measurement 
System  During  Experiments  on  Graphite  Specimens 


Graphite 

Temperature 

Specimen  Heating  Rate 

Radiative  Heat  Loss 

Grade 

K 

K's" 

i 

from  Specimen 
%  of  input  power 

"Fast" 

"Slow" 

"Fast" 

"Slow" 

1500 

6000 

4400 

2 

3 

Poco 

(Specimen^!) 

2000 

6000 

4000 

7 

10 

2500 

6900 

3800 

13 

22 

3000 

6200 

2700 

26 

45 

1500 

3800 

2900 

2 

3 

Py rolytic 

2000 

5100 

3600 

4 

6 

(Specimen-3) 

2500 

6600 

5200 

7 

9 

3000 

7500 

5700 

12 

15 
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The  heat  capacity  results  obtained  for  the  "fast"  and  "slow" 
experiments  were  fitted  separately  for  each  specimen  by  polynomial 
functions  of  temperature  using  the  least  squares  method.    The  best  fits 
were  quadratic  for  Poco  graphite  and  linear  for  pyrolytic  graphite, 
with  standard  deviations  (of  an  individual  point)  in  the  range  0.4  to 
0.67o.     A  study  of  the  results  indicated  the  following:     (a)  the 
difference  in  the  results  corresponding  to  "fast"  and  slow  experiments 
for  a  given  specimen  (in  the  range  0.2  to  0.6%)  is  less  than  or 
comparable  to  the  measurement  resolution,  and  (b)  the  difference  in 
the  results  corresponding  to  the  two  specimens  of  the  same  graphite  grade 
(0.37,  for  Poco  and  0.57,  for  pyrolytic)  is  less  than  or  comparable  to 
the  measurement  resolution.     It  may  be  concluded  that  for  a  given 
graphite  grade,   two  different  specimens  and  two  different  heating  rates 
yielded  the  same  value  for  heat  capacity.     However,    the  average  absolute 
difference  between  the  results  of  Poco  and  pyrolytic  graphites  was  1.87,; 
the  results  of  Poco  being  higher  than  those  of  pyrolytic. 

The  functions  for  heat  capacity  that  represent  the  combined  results 
of  the  two  specimens  and  the  two  heating  rates  for  each  graphite  grade  in 
the  range  1500  to  3000  K  are: 

For  Poco  graphite  (standard  deviation,  0.57,) 

Cp  =  19.84  +  3.490  x  10"3T-4.119  x  lO'V  (I) 

For  pyrolytic  graphite  (standard  deviation,  0.67.) 

C    =  21.60  +  1.534  x  10"3T  ) 
P 

where  T  is  in  K  and  C    is  in  J.mol"«K     .     In  the  computations  of  heat 

P 

capacity,   the  atomic  weight  of  graphite  was  taken  as  12.011.     The  heat 

8 


Table  3 

Heat  Capacity  of  Poco  and  Pyrolytic 
Graphites  According  to  Equations  (1)  and  (2) 


Temperature 
K 

Heat  Capacity,  J-mol  »K 

Poco 

Pyrolytic 

1500 

24.15 

23.90 

1600 

24.37 

24.05 

1700 

24.58 

24.21 

1800 

24.79 

24.36 

1900 

24.98 

24.51 

2000 

25.17 

24.67 

2100 

25.35 

24.82 

2200 

25.52 

24.97 

2300 

25.69 

25.13 

2400 

25.84 

25.28 

2500 

25.99 

25.44 

2600 

26.13 

25.59 

2700 

26.26 

25.74 

2800 

26.38 

25.90 

2900 

26.50 

26.05 

3000 

26.60 

26.20 
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capacity  for  the  two  grades  of  graphite  computed  using  the  Equations 

(1)  and  (2)  is  presented  in  Table  3.     Figures  1  and  2  show  the 
deviation  of  the  measured  heat  capacity  values  of  Poco  and  pyrolytic 
graphite  specimens  from  the  values  calculated  from  Equations  (1)  and 

(2)  ,  respectively. 

4.2.     Electrical  Resistivity 

Electrical  resistivity  of  the  graphite  specimens  was  determined 
from  the  same  experiments  that  were  used  to  calculate  heat  capacity. 
In  the  case  of  pyrolytic  graphite,   the  results  correspond  to  the 
measurements  in  the  direction  parallel  to  the  basal  planes.    The  differences 
in  the  results  between  "fast"  and  "slow"  experiments  for  a  given 
specimen  were  within  the  measurement  resolution.    However,  unlike  heat 
capacity,   there  were  differences  in  the  electrical  resistivity  of  the 
two  specimens  belonging  to  the  same  graphite  grade.     The  average  absolute 
difference  between  the  two  specimens  was  about  7%  for  Poco  graphite  and 
1.57o  for  pyrolytic  graphite.     It  should  be  noted,  however,   that  the  Poco 
graphite  specimens  represented  two  different  types,   therefore,  a 
difference  in  their  electrical  resistivity  may  be  expected.  The 
difference  in  the  case  of  pyrolytic  graphite  may  largely  be  attributed  to 
the  differences  in  the  fabrication  of  the  specimens.    The  electrical 
resistivity  of  Poco  graphite  at  temperatures  above  1500  K  was  about  four 
times  greater  than  that  of  pyrolytic  graphite. 
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The  functions  for  electrical  resistivity  that  represent  the 
average  of  the  results  corresponding  to  two  different  heating  rates  for 
the  two  types  of  Poco  graphite  specimens  in  the  range  1500  to  3000  K 
are : 

Poco  (AXM-5Q)  graphite  (standard  deviation,  0.07  %  ) 

p  =  908.56-1.4324  x  lO^T  +  1.4766  x  10_AT2-2 .2603  x  10_8T3  (3) 

Poco  (DFP-2)  graphite  (standard  deviation,  0.05%) 

D  =  808.97-1.0082  x  10-1T  +  1.2824  x  lO^I2 -1 . 9763  x  10"8T3  (4) 

The  function  for  electrical  resistivity  that  represents  the 
average  of  the  results  for  the  two  pyrolytic  graphite  specimens 
corresponding  to  two  different  heating  rates  in  the  range  1500  to 
3000  K  is: 

Pyrolytic  graphite  (standard  deviation,  0.9%) 

P  =  190.21-1.3574  x  10"3T  +  1.6492  x  lO'V  (5) 

-8 

where  T  is  in  K  and  p  is  in  10      Q«m.    The  electrical  resistivity  of 

the  graphites  computed  using  Equations  (3),   (4)  and  (5)  is  presented  in 

Table  4.    The  measured  electrical  resistivity  values  of  the  graphite 

specimens  are  shown  in  Figure  3.     The  data  on  both  Poco  graphite  specimens 

indicate  the  presence  of  an  inflection  point  in  the  function  of  electrical 

resistivity  versus  temperature  in  the  vicinity  of  2170  K  (2178  K  for 

specimen -1  and  2163  K  for  specimen-2) .    The  data  on  pyrolytic  graphite 

do  not  indicate  any  inflection  point  in  the  present  temperature  range. 

Electrical  resistance  measurements  at  283  K  with  a  Kelvin  bridge 

before  and  after  the  pulse  experiments  were  in  agreement  within  0.1%  for 

the  Poco  graphite  specimens  and  0.5%  for  the  pyrolytic  graphite  specimens. 

This  indicates  that  no  significant  irreversible  structural  or  chemical 
changes  had  taken  place  in  the  specimens  as  the  result  of 
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Table  4 

Electrical  Resistivity  of  Poco  and  Pyrolytic 
Graphites  According  to  Equations  (3),   (4)  and  (5) 


Temperature 
K 

-8 

Electrical  Resistivity,   10  Q-m 

Poco 

Pyrolytic 

Specimens 
3  and  4 

AXM-5Q 
(Specimen-1) 

DFP-2 
(Specimen-2) 

1 DUU 

949.7 

879.6 

225.3 

i  £aa 

964.8 

895.0 

230.3 

i  ir\f\ 
i  /UU 

980.7 

911.1 

235.6 

i  qaa 
loUU 

997.3 

927.7 

241.2 

1014.4 

944  .8 

247.2 

o  aaa 

1031.9 

962.2 

253.5 

1  1  AA 

2  100 

1049.6 

979.8 

260. 1 

1  1  A  A 

zzUO 

1067.4 

997.4 

267.0 

O  O  A  A 

1085.2 

1015.0 

274.3 

O/.  AA 

1  inn  o 

1 102  . 8 

1032 . 5 

281.9 

0  £  AA 

1120.2 

1049.6 

289.9 

2600 

1137.1 

1066.4 

298.2 

2700 

1153.4 

1082.6 

306.8 

2800 

1169.0 

1098.2 

315.7 

2900 

1183.7 

1113.1 

325.0 

3000 

1197.5 

1127.1 

334.6 
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Figure  3.        Electrical  resistivity  of  Poco  and  pyrolytic  graphites 
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their  having  been  heated  to  high  temperatures.     The  average  values  of  the 

-8 

electrical  resistivity  (in  10      Q'To)  for  the  graphites  at  283  K  were 
as  follows:     1475  for  Poco  (AXM-5Q),   1338  for  Poco  (DFP-2)  and  511 
for  pyro lytic  graphite. 

4.3.    Hemispherical  Total  Quittance 

Hemispherical  total  emittance  was  computed  for  each  specimen 
using  data  taken  during  both  heating  and  cooling  periods.     In  the 
case  of  pyrolytic  graphite,   the  results  correspond  to  radiation  emitted 
from  the  surface  parallel  to  the  basal  planes.     The  differences  between 
the  "fast"  and  "slow"  experiments  for  a  given  specimen  were  within 
the  measurement  resolution.    However,   the  results  were  different  for 
different  specimens.    The  functions  for  hemispherical  total  emittance 
that  represent  the  results  of  the  graphites  are: 

Poco  (AXM-5Q)  graphite,  1800  to  2900  K  (standard  deviation,  0.7%) 
e  =  0.679  +  6.00  x  10"5  T  (6) 
Poco  (DFP  -  2)  graphite, 1700  to  2900  K  (standard  deviation,  0.7%) 
e  =  0.794  +  2.28  x  10"5  T  (7) 
Pyrolytic  graphite,    2300  to  3000  K  (standard  deviation,  1.2%) 
e  =  0.641-5.70  x  10"5  T  (8) 
where  T  is  in  K.     The  hemispherical  total  emittance  of  the  graphites 
computed  using  Equations  (6),   (7)  and  (8)  is  presented  in  Table  5.  The 
measured  values  of  the  hemispherical  total  emittance  of  the  graphite 
specimens  are  shown  in  Figure  4. 
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Table  5 

Hemispherical  Total  Emittance  of  Poco 
and  Pyrolytic  Graphites  According  to 


Equations 

(6), 

(7) 

and  (8) 

Hemispherical  Total  Emittance 

Temperature 

Poco 

Pyrolytic 

K 

AXM-5Q 

DFP-  2 

(Specimen 

-1) 

(Specimen-2) 

(Specimen-3) 

1700 

0.833 

1800 

0.787 

0.835 

1900 

0.793 

0.837 

2000 

0.799 

0.840 

2100 

0.805 

0.842 

2200 

0.811 

0.844 

2300 

0.817 

0.846 

U  .  J  IV 

2400 

0.823 

0.849 

U  .  JUH 

2500 

0.829 

0.851 

0.499 

2600 

0.835 

0.853 

0.493 

2700 

0.841 

0.856 

0.487 

2800 

0.847 

0.858 

0.481 

2900 

0.853 

0.860 

0.476 

3000 

0.470 
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5.     Estimate  of  Errors 

Estimates  of  errors  in  measured  and  computed  quantities  lead  to 
the  following  estimates  of  errors  in  the  properties:     heat  capacity, 
3%;  electrical  resistivity,   1%  (Poco)  and  3%  (pyrolytic) ;  hemispherical 
total  emittance,  5%.     Details  regarding  the  estimates  of  errors  and 
their  combination  in  experiments  using  the  present  measurement  system 
are  given  in  a  previous  publication  [4].     Specific  items  in  the  error 
analyses  were  recomputed  whenever  the  present  conditions  differed  from 
those  in  the  earlier  publication. 

Because  of  the  high  value  of  the  hemispherical  total  emittance 
of  graphite  (especially  Poco  grades),  heat  loss  from  the  specimen  due 
to  thermal  radiation  constitutes  a  substantial  fraction  of  the  input 
power,  especially  at  temperatures  above  2500  K  (Table  2).  However, 
as  a  result  of  data  taken  during  the  cooling  period,  an  accurate  account 
of  its  magnitude  in  correcting  heat  capacity  was  possible.     This  is 
demonstrated  by  the  rather  closely  agreeing  results  (within  measure- 
ment resolution)  of  the  two  sets  of  experiments  performed  on  the  same 
specimen  at  two  different  heating  rates. 

6.  Discussion 

6.1.    Heat  Capacity 
Only  a  limited  number  of  experimental  heat  capacity  results  for 
graphite  at  temperatures  above  1500  K  was    found  in  the  literature. 
Since  most  of  them  are  for  different  grades  of  graphite,  it  is  question- 
able whether  one  has  a  common  ground  for  comparison.    As  it  may  be  seen 

in  Figure  5,  considerable  differences  exist  between  the  results  of 
various  investigators.  ,Q 


1000  1500  2000  2500  3000 

TEMPERATURE,  K 

Figure  5.        Heat  capacity  of  graphite  reported  in  the  literature. 

The  measurements,  with  the  exception  of  those  of  Rasor 
and  McClelland  and  the  present  work,  were  performed 
using  drop  calorimetric  techniques. 


The  heat  capacity  results  of  the  present  work  agree    best  with 
those  reported  by  West  and  Ishihara  [6]  on  grade  CCH  graphite.    In  the 
overlapping  temperature  region  (1500  to  2600  K),  the  results  of  West 
and  Ishihara  are,  on  the  average,  about  0.6%  lower  than  the  Poco 
graphite  results  and  about  1.2%  higher  than  the  pyro lytic  graphite 
results.    This  rather  close  agreement  is  even  more  significant  if 
one  considers  that  different  measurement  techniques  (drop  and  pulse) 
were  employed  in  the  two  cases. 

The  smoothed  results  of  Rasor  and  McClelland  [7]  are  average 
values,  as  given  by  the  authors,   representing  data  on  four  grades  of 
graphite  ^3474D,  7087,  GBH,  and  GBE) .    Their  actual  data  have  a  scatter 
of  approximately  5%  about  their  smooth  curve.    The  difference  between 
the  results  of  Rasor  and  McClelland  and  those  of  the  present  work  is, 
on  the  average,  about  3%  (Poco)  and  5%  (pyrolytic). 

Lucks  et  al.   [8]  have  reported  results  for  two  grades  of  graphite. 
The  upper  and  lower  curves  (in  Figure  5)  represent  the  results  for  the 
grades  7087  and  GBH,   respectively.    The  disagreement  between  the  results 
of  Lucks  et  al.  and  Rasor  and  McClelland  for  the  same  grades  of  graphite 
is  more  than  10%,  suggesting  serious  experimental  errors. 

In  the  temperature  region  of  overlap  with  the  present  work  (1500  to 
1800  K),  the  results  for  ATJ  graphite  reported  by  Fieldhouse  et  al.  [9] 
agree,  within  3%,  with  those  of  the  present  work.    However,   the  slope  of 
the  curve  representing  the  results  of  Fieldhouse  et  al.  is  considerably 
greater  than  the  slopes  of  the  curves  representing  the  present  work. 
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The  results  reported  recently  by   Sheindlin  et  al.  [10]  on  a 
number  of  different  graphites  fall  between  the  present  work  results  on 
the  two  graphite  grades  in  the  range  1500  to  2300  K.    However,  above 
2300  K,  their  results  indicate  a  sharp  increase  in  heat  capacity 
which  reaches  a  value  of  29.55  J'mol"1^"1  at  3000  K,  which  is  11%  higher 
than  the  present  work  result  on  Poco  graphite. 

In  addition  to  the  present  measurements,   the  only  other  known 
measurements  on  Poco  (AXM-5Q)  graphite  were  those  reported  in  the  AGARD 
Report  [11]  covering  the  temperature  range  373  to  1273  K.    Since  these 
data  have  no  overlapping  region  with  the  present  results  it  is  difficult 
to  estimate  the  differences  in  the  two  results.    However,   from  Figure 
5  it  may  be  seen  that  an  extrapolation  of  the  literature  results  to 
1500  K  is  likely  to  yield  a  value  approximately  2%  lower  than  that  of 
the  present  work  on  Poco  graphite. 

The  rather  closely  agreeing  results  of  West  and  Ishihara  [6]  and 
those  of  the  present  work  were  obtained  on  different  grades  of  graphite 
using  different  measurement  techniques  indicating  that  the  heat 
capacity  of  graphite  may  not  be  very  sensitive  to  differences  in  certain 
grades.    Therefore,   the  differences  in  the  literature  results  should  not 
be  entirely  attributed  to  grade  differences,  and  the  foregoing  observation 
suggests  that  these  differences  might  have  been  partly  due  to  experimental 
errors.     In  addition,   the  state  of  the  specimen  (purity,  degree  of 
graphitization  and  annealing)  have  undoubtedly  played  an  important  role 
in  the  results.     Sheindlin  et  al.  [10]  have  also  arrived  to  a  similar 
conclusion  as  a  result  of  their  measurements  on  several  graphites.  They 
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have  observed  that  the  values  obtained  for  different  graphites  did 
not  differ  by  more  than  their  experimental  error  (1.5  to  2.5%  in 
enthalpy)  . 

It  may  be  seen  that  heat  capacity  of  graphite  reaches  the  Dulong 
and  Petit  value  of  3R  at  a  temperature  in  the  vicinity  of  2000  K,  and 
continues  to  increase  above  that  temperature. 

The  next  challange  in  heat  capacity  measurements  is  likely  to  be 
at  temperatures  above  3000  K  where  measurements  will  yield  information 
of  practical  and  theoretical  significance.     It  is  important  to  resolve  the 
questions  of  whether  the  heat  capacity  of  graphite  reaches  a  constant 
value  at  high  temperatures,  or  increases  sharply,  as  was  reported  by 
Rasor  and  McClelland  [7]  and  Sheindlin  et  aL  [10]. 

6.2.     Electrical  Resistivity 

A  comparison  of  the  present  work  and  literature  results  on  the 
electrical  resistivity  of  the  two  graphite  grades  is  presented  in 
Figure  6  . 

In  addition  to  the  present  measurements,   the  only  other  known 
measurements  of  the  electrical  resistivity  of  Poco  graphite  above 
1500  K  are  those  reported  in  the  AGARD  Report  [11].    A  comparison  of  the 
present  work  results  with  the  two  different  measurements  indicates  that 
in  one  case  the  present  work  values  are  8%  lower  (at  1500  K)  and  1%  higher 
(at  2500  K)  than  the  literature  values;  in  the  other  case,  an  extra- 
polation over  200  K  indicates  that  at  1500  K  the  present  work  value  is 
about  5%  higher  than  the  literature  value. 
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Figure  6.        Electrical  resistivity  of  Poco  (AXM-5Q)  and  pyro lytic 
graphite  reported  in  the  literature.     The  AGARD  Work" 


corresponds  to  measurements  by. di-^e rent  participants 
reported  in  Reference  [11]. 
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Electrical  resistivity  of  pyrolytic  graphite  was  reported  by 
Pappis  and  Blum  [121  up  to  1700  K  and  by  Petrov  et  al.   [13]  up  to 
2200  K.     At  1500  K  the  present  work  result  is  2  and  3%  lower  than  the 
values  reported  by  Pappis  and  Blum,  and  Petrov  et  al.,  respectively; 
at  1700  K  it  is  TL  lower  than  that  of  Pappis  and  Blum;  at  2200  K  it 
is  5%  higher  than  that  of  Petrov  et  al. 

6.3.     Hemispherical  Total  Emittance 

No  other  measurements  of  the  hemispherical  total  emittance  of 
Poco  graphite  have  been  located  in  the  literature.     The  results  for 
other  graphites  such  as  Acheson  graphite  reported  by  Jain  and  Krishnan 
[14]  yield  0.81  at  1700  K  and  0.83  at  2100  K  which  are  within  370  of 
the  results  of  the  present  work. 

Because  of  differences  in  the  methods  of  fabrication  for  the  specimens 
and  the  anisotropic  properties  of  pyrolytic  graphite,  considerable 
differences  in  hemispherical  total  emittance  measurements  by  different 
investigators  may  be  expected.    As  presented  by  Petrov  et  al.    [15],  the 
values  for  emittance  reported  in  the  literature  at  temperatures  above 
1200  K  range  from  0.5  to  0.9.     In  addition  to  the  large  differences, 
the  trends  (emittance  versus  temperature)  are  also  considerably  dif- 
ferent, in  some  cases  emittance  is  reported  to  increase  with  increasing 
temperature  and  in  others  to  decrease  with  temperature.     The  results 
of  the  present  work  indicate  a  decrease  of  emittance  with  temperature 
and  their  closest  agreement  is  with  those  of  Petrov  et  al .    [15].  The 
measurements  by  Petrov  et  al.  show  a  decrease  in  emittance  up  to  about 
2000  K  and  a  small  increase  above  2000  K.     Extrapolation  of  the  present 
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work  results  to  2000  K  yields  a  value  of  0.53  which  is  about  5%  lower 
than  the  value  given  by  Petrov  et  al. 
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Chapter  2 


SIMULTANEOUS  MEASUREMENTS  OF  HEAT  CAPACITY, 
ELECTRICAL  RESISTIVITY  AND  HEMISPHERICAL  TOTAL 
EMITTANCE  BY  A  PULSE  HEATING  TECHNIQUE: 
VANADIUM,   1500  TO  2100  K* 

A.  Cezairliyan,  F.  Righini**  and  J.  L.  McClure 
Institute  for  Materials  Research 
National  Bureau  of  Standards 
Washington,  D.  C.  20234 

Abstract 

Simultaneous  measurements  of  heat  capacity,  electrical 
resistivity  and  hemispherical  total  emittance  of  vanadium  in  the 
temperature  range  1500  to  2100  K  by  a  subsecond  duration,  pulse 
heating  technique  are  described.     The  results  are  expressed  by  the 

relations : 

c    =  56.34  -  3.839  x  10~2T  +  1.563  x  10_5T2 
P 

p  =  8.794  +  6.282  x  10_2T  -6.804  x  10_6T2 

- 1  - 1  -8 
where  c^  is  in    j-mol    -K    ,  p  is  in  10     Q»m,  and  T  is  in  K.  The 

values  for  the  hemispherical  total  emittance  are:     0.313  at  1900  K 

and  0.332  at  2000  K.     Estimated  inaccuracies  of  the  measured  propert 

are:     3%  for  heat  capacity,  0.5%  for  electrical  resistivity  and 

5%  for  hemispherical  total  emittance. 
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1.     Introduction  , 

In  this  paper,  application  of  a  pulse  heating  technique  to 
the  simultaneous  measurements  of  heat  capacity,  electrical  resistivity 
and  hemispherical  total  emittance  of  vanadium  in  the  temperature 
range  1500  to  2100  K  is  described. 

The  method  is  based  on  rapid  resistive  self-heating  of  the 
specimen  from  room  temperature  to  high  temperatures  (above  1500  K) 
in  less  than  one  second  by  the  passage  of  an  electrical  current 
pulse  through  it;  and  on  measuring,  with  millisecond  resolution, 
such  experimental  quantities  as  current  through  the  specimen,  potential 
drop  across  the  specimen,  and  specimen  temperature.     Details  regarding 
the  construction  and  operation  of  the  measurement  system,   the  methods 
of  measuring  experimental  quantities,  and  other  pertinent  information, 
such  as  the  formulation  of  relations  for  properties,  error  analysis, 
etc.  are  given  in  earlier  publications  [1,  2]^". 

In  the  following  sections  of  this  paper  a  new  approach  -  tabular 
format  -  is  adopted  in  presenting  information  on  the  specimen,  measure- 
ments, system  characteristics,  results,  and  errors. 

The  reasons  for  adopting  this  format  are:     (1)  to  facilitate  the 
preparation  of  manuscripts,   (2)  to  standardize  the  contents  of  the 
papers  on  measurements  using  the  present  or  similar  systems  in  other 
laboratories,   (3)  to  facilitate  information  retrieval  by  the  reader, 
(4)  to  provide  a  means  for  efficient  and  accurate  identification, 


Figures  in  brackets  indicate  the  literature  references  at  the  end 
of  this  paper. 
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TABLE  1 


Specimen  information 


No. 

Item 

Unit 

Explanation 

i 

i 

Subs  tance 

: 

Vr  a  0  da  t  urn  ^po  iy  cryb  Lai  line  ) 

Z 

* 

Source 

Materials  ttesearcn  corporation 

3 

Purity 

99 . 9  % 

4 

Impurities 

Listed  in  table  la 

5 

Geomet  ry 

Tube  made  from  rod  by  electro-erosion 

6 

Dimensions 

total  length 

mm 

/O  .  ZD 

•k-k 

effective  length 

mm 

25.53 

outside  diameter 

mm 

6.3 

wall  thickness 

mm 

0.5 

blackbody  hole 

mm 

0.5  x  1 

7 

Weight 

total  weight 

e 

o 

4.319 

"kit 

effective  weight 

g 

1.438 

8 

Characteristics 

^  1"  nm i  p  up  i  oVi  t- 

density 

-3 

g-cm 

6.1 

resistivity  at  293  K 

10"8Q-m 

21.72 

9 

Special  treatment 

Heat  treated  by  pulse  heating  before 

the  experiments  -  30  pulses  to  1900  K 

* 

The  supplier  is  identified  in  this  paper  in  order  to  adequately  characterize  the 
specimen  -  Such  an  identification  does  not  imply  recommendation  or  endorsement 
by  the  National  Bureau  of  Standards. 

** 

Effective  refers  to  the  portion  of  the  specimen  between  the  voltage  probes. 
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coding,  characterization  and  data  reduction  by  scientific  and 

J 

technical  information  centers,  and  (5)  to  ultimately  lead  to  the 
computerized  preparation  of  papers. 

2.  Measurements 

The  details  regarding  the  vanadium  specimen  used  in  the  present 
measurements  are  given  in  table  1.    A  summary  of  the  measurement 
technique  and  the  operational  characteristics  of  the  system  is  given 
in  table  2.     The  polynomial  functions  (obtained  by  the  least  squares 
method)  that  represent  the  experimental  results  are  given  in  table  3. 
The  final  values  on  properties  at  100  degree  temperature  intervals 
computed  using  the  functions  are  given  in  table  4.    The  experimental 
results  are  presented  in  the  Appendix.     Each  number  tabulated  in 
the  Appendix  represents  results  from  over  fifty  original  data  points. 
An  estimate  of  errors  in  the  measured  and  computed  quantities  is  given 
in  table  5.     All      values      reported  in  this  paper  are  based  on  the 
International  Practical  Temperature  Scale  of  1968  [3].     In  all 
computations,   the  geometrical  quantities  are  based  on  their  room 
temperature  (298  K)  dimensions. 


Table  la 


Element 

C 

Fe 

Nb 

N 

0 

P 

ppm 

120 

20 

60 

10 

15 

15 

Element 

Si 

Ta 

Ti 

W 

Zr 

ppm 

50 

70 

10 

30 

15 

a 

The  total  amount  of  all  other  detected  elements  is  less  than 
50  ppm,   each  element  being  below  10  ppm  limit. 
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TABLE  2 


Measurement  technique  and  system  characteristics 


No. 


I  tern 


Unit 


Explanation  and  Data 


1  General  technique 

2  ;     Voltage  measurement 

3  Current  measurement 

i 

4  j  Temperature  measurement 

5  Specimen  environment 

6  Power  source 

Recording 
Signal  resolution 
9  Time  resolution 

Data  processing 
Date  of  measurements 

12  Number  of  experiments 

13  J     Temperature  range 

14  |     Temperature  subranges 


Experiment  duration 

16  Current  pulse  length 

17  Imparted  power 

18  Current 

I 

19  i  Rate  of  current  change 

20  |  Heating  rate 

21  Cooling  rate 

22  Radiative  heat  loss 
 1        (7o  of  input  power) 


ms 


K 
K 


ms 
ms 
W 
A 


A-ms 


1 


K-ms 


K-ms 


-1 
-1 


Pulse  heating  (subsecond)    [1,  2] 

Across  tungsten  knife-edge  probes 

Across  standard  resistor  (0.001  Q) 
in  series  with  the  specimen 

High-speed  photoelectric  pyrometer  [4] 


Vacuum 


1.3  x  10"3  N-m"2  (- 


10  torr) 


Battery  bank  (14  series-connected 
2V  batteries,   capacity  1100  Ah  each) 

Digital  data  acquisition  system 

~  0.01%  (at  full  scale) 

0.4 

Time -sharing  computer 

August  1973 

4 

1500  -  2100 

I  (1480  -  1710) 

II  (1690  -  1890) 

III  (1810  -  2050) 

IV  (2000  -  2130) 

600  -  620 

400  -  420 

4000  -  5200 

1300  -  1600 

0.3  -  0.6 

3.9  -  4.5 

0.04  -  0.15 

1%  at  1500  K 
47o  at  2100  K 
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TABLE  3 


Functional  representation  of  results  on    va-O&ci  Lum 


Heat  capacity   (J-mol  "*"K  _) 


-8, 


Resistivity  (10  Q-m) 


Hemispherical  total  emittance 


c    =  A  +  BT  +  CT 

P 

\  =  56.34 


B  =  -3.839  x  10 
C  =  1.563  x  10"5 
1500  K  <  T  <  2100  K 
*  =  1.1% 


o  =  A  +  BT  +  CT 
A  =  8.794 
B  =  6.282  x  10"2 
C  =  -6.804  x  10"6 
1500  K  <  T  <  2100  K 
CT    =  0.06% 


e  =  A  +  BT 

A  =  -5.413  x  10"' 

B  =  1.930  x  10 


1880  K  <  T  <  2050  K 
a*  =  0.05% 


Standard  deviation  of  an  "individual  point' 


TABLE  4 

Results  on  properties  of   N/3,  na  ol  iurn 


T 
(K) 

cp 

(J -mol 

p 

(10~* 

e 

1500 

33 

.92 

87 

72 

1600 

34 

.93 

91 

89 

1700 

36 

.25 

95 

92 

1800 

37 

.88 

99 

83 

1900 

39 

.82 

103 

59 

0.313 

2000 

42 

08 

107 

22 

0.332 

2100 

44 

.65 

110 

71 

35 


TABLE  5 


Error  analysis  (at  2000  K) 


Quantity 

★ 

Imprecision 

** 

Inaccuracy 

Temperature 

0.5  K 

4  K 

Voltage 

0 . 037. 

0.17. 

Current 

u .  UJ7. 

U.J./. 

Heat  capacity 

17. 

37. 

Electrical  resistivity 

0.17. 

0.57. 

Hemispherical  total 

emi  ttance 

0.17. 

57. 

Imprecision  refers  to  the  standard  deviation  of  an  individual 
point  as  computed  from  the  difference  between  measured  value  and 
that  from  the  smooth  function  obtained  by  the  least  squares  method. 

* 

Inaccuracy  refers  to  the  estimated  total  error  (random  and  systematic). 
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3.  Discussion 

The  heat  capacity,   electrical  resistivity  and  hemispherical  total 
emittance  of  vanadium  measured  in  this  work  are  presented  and  compared 
graphically  with  those  reported  in  the  literature  in  figures  1,   2  and 
3,  respectively. 

The  heat  capacity  results  of  this  work  are  approximately  6-8% 
lower  than  those  of  Jaeger  and  Veenstra  [5]  and  Fieldhouse  and  Lang  [6], 
and  are  approximately  8%  higher  than  those  of  Peletskii  et  al.  [7] 
over  the  respective  overlapping  regions.     The  results  reported  in  the 
literature  were  for  temperatures  below  1900  K.     In  this  work,  the 
measurements  are  extended  to  2100  K,  which  is  approximately  90  K  below 
the  melting  point  of  vanadium. 

The  electrical  resistivity  results  are  in  reasonably  good  agreement 

(maximum  differences  about  270)  with  those  reported  by  Horz  et  al.  [8] 

and  Peletskii  et  al.   [71.     At  293  K,   the  electrical  resistivity  of  this 

-8 

work  is  approximately  37o  higher  than  the  value  (21.02  x  10  Q-m) 
reported  by  Peletskii  et  al.   [7].     In  the  range  of  the  present  measurements, 
vanadium  showed  a  negative  departure  from  linearity  in  the  temperature 
dependence  of  electrical  resistivity.     A  similar  behavior  was  previously 
noted  for  niobium  [9]  and  tantalum  [10],  which  belong  to  the  same  group 
(V)  on  the  periodic  table. 

Considerable  differences  in  hemispherical  total  emittance  results 
of  various  investigators  may  be  expected  due  to  the  differences  in 
specimen  surface  conditions.     The  difference  between  the  present  results 
and    the  extrapolated  values  of  Horz  et  al.   [8]  is  approximately  6%. 
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FIGURE  1.        Heat  capacity  of  vanadium  as  reported  in  the  literature, 
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As  was  the  case  with  the  earlier  results  on  other  refractory  metals  [2, 

9,   10,  11]  obtained  with  the  present  measurement  system,  heat  capacity 

of  vanadium  at  high  temperatures  is  considerably  higher  than  the 

Dulong  and  Petit  value  of  3R.     Some  of  this  departure  is  due  to 

c  -c    and  the  electronic  terms.    However,   they  do  not  account  for 
P  v 

the  entire  difference.    Heat  capacity  above  the  Debye  temperature 
may  be  expressed  by 

c„  =  A-^+CT+bc 

-1  -1  -2 

where  the  constant  term  is  3R  (24.943  J.mol  »K    ),  the  term  in  T 

is  the  first  term  in  the  expansion  of  the  Debye  function,   the  term 
in  T  represents  Cp~cv  and  electronic  contributions,  and  the  quantity 
Ac  represents  excess  in  measured  heat  capacity  at  high  temperatures, 

which  is  not  accounted  for  by  the  first  three  terms.    The  coefficients 

5  -3 
B  (1.50  x  10  )  and  C  (4.86  x  10    )  were  obtained  from  heat  capacity 

values  at  270  K  [12,   13]  and  at  1000  K.    The  value  at  1000  K  was 

estimated  based  on  the  extrapolation  of  the  present  work  results 

and  the  low  temperature  results  reported  in  the  literature  [12,  13]. 

Using  the  above  equation  and  the  heat  capacity  results  of  this 
work,  the  quantity  Ac  was  computed  for  temperatures  above  1500  K. 
The  results  for  Ac  in  J.mol'^K-1  are:     1.8  at  1500  K,  4.2  at  1800  K, 
and  9.5  at  2100  K. 
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Although  the  mechanisms  of  vacancy  generation  become  important 
at  high  temperatures,  it  was  not  possible  to  attribute  the  high  heat 
capacity  values  entirely  to  vacancies.     To  demonstrate  this,  a  crude 
estimate  of  the  contribution  of  vacancies  to  heat  capacity  was  made 
using  the  method  described  in  a  previous  publication  [2],    The  results 
indicate  that  vacancy  contribution  would  be  small,  approximately 
0.03  J-mol'^K-1  at  1500  K  and  0.5  J-mol~»K~1  at  2100  K,  and  would  not 
account  for  the  high  heat  capacity  values.    If  the  entire  difference 
between  measured  and  computed  (using  the  first  three  terms  in  the 
above  equation)  heat  capacities  is  attributed  to  vacancies,  values  of 
1.1  eV  for  vacancy  formation  energy  and  3.87o  for  vacancy  concentration 
at  the  melting  point  are  obtained.     Both  of  these  values,  especially 
the  concentration,  seem  to  be  unrealistic  for  vanadium. 

The  new  format  adopted  in  this  paper  has  demonstrated  the 

feasibility  of  concise  yet  complete  and  systematic  documentation 
of  measurements  performed  with  the  present  system.     For  reasons 
enumerated  in  the  Introduction,   it  is  hoped  that  the  approach  of 
a  standardized  quasi -tabular  presentation  will  find  wide  acceptance 
among  the  experimenters  in  the  same  and  related  fields. 
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4.  Appendix 


TABLE  A-l 


Experimental  results  on  heat  capacity  and  electrical  resistivity  of  vanadium 


Temperature 
(K) 

(J  mol^K  ) 

* 

ACp 
(%) 

p 

(10"8n-m) 

* 

AP 
(%) 

1500 

33.48 

-  1.32 

87.66 

0.06 

1550 

34.50 

+  0.33 

89.81 

0.01 

1600 

35.32 

+  1.10 

91.93 

+ 

0.05 

1650 

35.88  > 

+  0.92 

94.03 

+ 

0.12 

1700 

35.81 

-  1.23 

95.83 

- 

0.09 

1750 

37.20 

+  0.46 

97.86 

- 

0.03 

1800 

38.26 

+  0.99 

99.87 

+ 

0.05 

1850 

38.61 

-  0.53 

101.69 

0.03 

1900 

39.71 

-  0.30 

103.56 

0.02 

1950 

40.80 

-  0.29 

105.40 

0.02 

2000 

41.96 

-  0.30 

107.20 

0.01 

2050 

42.72 

-  1.43 

109.01 

+ 

0.03 

2100 

45.34 

+  1.51 

110.70 

0.00 

The  quantities  ACp  and  Ap  are  percentage  deviations  of  the  individual 
results  from  the  smooth  functions  represented  by  the  pertinent  equations 
in  table  3. 
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TABLE  A-2 


Experimental  results  on  hemispherical  total  emittance  of  vanadium 


Temperature 
(K) 

e 

* 

Ae 

a) 

Temperature 
(K) 

e 

* 

Ae 
(%) 

1874.5 

0.307 

0. 12 

2022.6 

0.336 

0.06 

1878.4 

0.308 

0.03 

2028.0 

0.337 

0.02 

1882 .4 

U .  JU9 

+ 

0.03 

ft  ^  i  a 

U  .  J  JO 

+ 

0.01 

1886.4 

0.310 

+ 

0.06 

2038.9 

0.339 

+ 

0.02 

1890.5 

0.311 

+ 

0.06 

2044.5 

0.341 

+ 

0.02 

1894.7 

0.312 

+ 

0.03 

2050.1 

0.342 

+ 

0.01 

*The  quantity  Ae  is  percentage  deviation  of  the  individual  results  from 
the  smooth  function  represented  by  the  pertinent  equation  in  table  3. 
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Chapter  3 

SIMULTANEOUS  MEASUREMENTS  OF  HEAT  CAPACITY, 
ELECTRICAL  RESISTIVITY  AND  HEMISPHERICAL  TOTAL 
EMITTANCE  BY  A  PULSE  HEATING  TECHNIQUE: 
ZIRCONIUM,   1500  TO  2100  K* 

A.  Cezairllyan  and  F.  Righini 
Institute  for  Materials  Research 
National  Bureau  of  Standards 
Washington,  D.  C.  20234 

Abstract 

Simultaneous  measurements  of  heat  capacity,  electrical 

resistivity  and  hemispherical  total  emittance  of  zirconium  in 

the  temperature  range  1500  to  2100  K  by  a  subsecond  duration,  pulse 

heating  technique  are  described.     The  results  are  expressed  by 

the  relations: 

C    =  36.65  -  1.435  x  10-2  T  +  6.624  x  10"6T2 
P 

o    =  87.95  +  1.946  x  10"2T 

e    =  0.2031  +  6.362  x  10"5T 

where  Cp  is  in  J.mol"1-  K~*     ,   o  is  in  10"8n -m,  and  T  is  in  K. 
Estimated  inaccuracies  of  the  measured  properties  are:     3%  for  heat 
capacity,  2%  for  electrical  resistivity  and  5%  for  hemispherical  total 
emi  ttance . 

This  work  was  supported  in  part  by  the  U.  S.  Air  Force  Office 
of  Scientific  Research. 

Guest  scientist  from  the  Istituto  di  Metrologia  "G.  Colonnetti" 
in  Torino,  under  a  fellowship  from  the  Consiglio  Nazionale  delle 
Ricerche  of  Italy. 


47 


I 

L.  Introduction 

In  this  paper,   application  of  a  pulse  heating  technique  to 
the  simultaneous  measurements  of  heat  capacity,  electrical 
resistivity  and  hemispherical  total  emittance  of  zirconium 
in  the  temperature  range  1500  to  2100  K  is  described. 

The  method  is  based  on  rapid  resistive  self-heating  of 
the  specimen  from  room  temperature  to  high  temperatures  (above 
1500  K)  in  less  than  one  second  by  the  passage  of  an  electrical 
current  pulse  through  it;  and  on  measuring,  with  millisecond 
resolution,   such  experimental  quantities  as  current  through  the 
specimen,  potential  drop  across  the  specimen,  and  specimen  temp- 
erature.    Details  regarding  the  construction  and  operation  of  the 
measurement  system,   the  methods  of  measuring  experimental  quantities, 
and  other  pertinent  information,   such  as  the  formulation  of  relations 
for  properties,   error  analysis,  etc.  are  given  in  earlier 
publications  [ 1,  2]^. 

In  the  following  sections  of  this  paper  a  tabular  format  is 
adopted  in  presenting  information  on  the  specimen,  measurements, 
system  characteristics,   results     and  errors. 


''"Figures  in  brackets  indicate  the  literature  references  at  the 
end  of  this  paper. 
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i  TABLE  1 
Specimen  information 


No. 

I  tern 

Unit 

Explanation 

1 

Substance 

Zirconium 

( po ly cry s  ta 1 1 ine ) 

2 

* 

Source 

Materials 

Research  Corporation 

3 

Purity 

99 

.98% 

4 

Impurities 

Listed  in 

table  la 

5 

Geometry 

Tube  made 

from  rod  by  electro-erosion 

6 

Dimensions  (nominal) 

total  length 

mm 

76 

.2 

-k-k 

effective  length 

mm 

25 

.4 

outside  diameter 

mm 

6 

.3 

wall  thickness 

mm 

0 

.5 

blackbody  hole 

mm 

0 

.5  x  1 

(rectangular) 

7 

Weight 

total  weight 

g 

4 

.312 

effective  weight 

g 

1 

.425 

8 

Characteris  tics 

atomic  weight     [ 4 1 

91 

.22 

■k-k-k 

density 

g.cm"3 

6 

.53 

resistivity  at  293  K 

10_8Q-m 

42 

.8 

The  supplier  is  identified  in  this  paper  in  order  to  adequately  characterize  the 
specimen  -  Such  an  identification  does  not  imply  recommendation  or  endorsement 
by  the  National  Bureau  of  Standards. 

■k-k 

Effective  refers  to  the  portion  of  the  specimen  between  the  voltage  probes. 

■k-k-k 

Value  measured  in  the  present  work. 
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Measurements 


The  details  regarding  the  zirconium  specimens  used  in  the 
present  measurements  are  given  in  table  1.    A  summary  of  the 
measurement  technique  and  the  operational  characteristics  of  the 
system  is  given  in  table  2.     The  polynomial  functions  (obtained  by 
the  least  squares  method)  that  represent  the  experimental  results 
are  given  in  table  3.     The  values  of  properties  at  100  degree 

temperature  intervals  computed  using  the  functions  are  given  in 
table  4.    The  experimental  results  are  presented  in  the  Appendix. 
Each  number  tabulated  in  the  Appendix  represents  results  from  over 
fifty  original  data  points.     An  estimate  of  errors  in  the  measured 
and  computed  quantities  is  given  in  table  5.    All  values  reported 
in  this  paper  are  based  on  the  International  Practical  Temperature 
Scale  of  1968  [3] .     In  all  computations,   the  geometrical  quantities 
are  based  on  their  room  temperature  (298  K)  dimensions. 


Table  la 
Impurities  in  the  specimen 


(according  to  the  manufacturer's  analysis) 


Element 

c 

H 

I 

0     i  N 

Al 

Fe 

Hf 

Ni 

Si 

Ti 

ppm 

6 



3.3 

125 

r  1 
2.1 

3 

30 

40 

1.5 

1.5 

1 

V 

The  total  amount  of  all  other  detected  elements  is  less  than  6  ppm, 
each  element  being  below  1  ppm  limit. 
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TABLE  2 


Measurement  technique  and  system  characteristics 


No. 

I  tern 

Unit 

Explanation  and  Data 

1 

General  technique 

Pulse  heating  (subsecond) 

2 

Voltage  measurement 

Across  tungsten  knife-edge  probes 

3 

Current  measurement 

Across  standard  resistor  (0.001  Cl) 
in  series  with  the  specimen 

4 
5 

Temperature  measurement 
Specimen  environment 

High-speed  photoelectric  pyrometer[5] 
Vacuum  ~1. 3  x  10  J  N-m  (~10_:>torr) 

6 

Power  source 

Battery  bank 
2V  batteries 

(14  series -connected 
,   capacity  1100  A*h  each) 

7 

Recording 

Digital  data 

acquisition  system 

8 

Signal  resolution 

~  0. 

01%  (at 

full  scale) 

9 

Time  resolution 

ms 

0.4 

10 

Data  processing 

Time 

-sharing 

computer 

11 

Number  of  specimens 

3 

12 

Number  of  experiments 

12 

13 

Temperature  range 

K 

1500 

-  2100 

14 

Temperature  subranges 

K 

I 
II 
III 
IV 

(1450  - 
(1680  - 
(1810  - 
(1840  - 

1680) 
1900) 
2050) 
2110) 

15 

Experiment  duration 

ms 

550 

-  680 

16 

Current  pulse  length 

ms 

350 

-  480 

17 

Imparted  power 

W 

1700 

-  3100 

18 

Current 

A 

730 

-  840 

19 

Rate  of  current  change 

A 

-1 

•  ms 

0.08 

-  0.13 

20 

Heating  rate 

K 

-1 

•ms 

-1 

•  ms 

2.9 

-  4.8 

21 

Cooling  rate 

K 

0.08 

-  0.3 

22 

Radiative  heat  loss 
(%  of  input  power) 

2%  at 
10%  at 

1500  K 
2100  K 

If 


I 
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TABLE  3 

Functional 

representation  of  results  on 

zirconium 

Heat  capacity  (J-mol  .^K  ^) 

-8 

Resistivity  (10  Q-m) 

Hemispherical  total  emittanrp 

2 

C    =  a  +  bT  +  cT 
P 

p  =  a  +  bT 

e  =  a  +  bT 

a  =  36.65 

a  =  87.95 

a  =  0.2031 

b  =  -1.435  x  10"2 

b  =  1.946  x  10"2 

b  =  6.362  x  10"5 

c  =  6.624  x  10"6 

1500  K  <  T  <  2100  K 

1500  K  <  T  <  2100  K 

1650  K  <  T  <  2050  K 

a*  =  0.7% 

a*  =  0.6% 

a*  =  0.9% 

Standard  deviation  as  computed  from  the  difference  between  the  value  of  an 
experimental  result  (as  tabulated  in  the  Appendix)  and  that  from  the  smooth 
functions  reported  above. 

TABLE  4 


Results  on  properties  of  zirconium 


T 

cp 

P 

e 

(K) 

(J-mol"1. K_1) 

(10"8Q.m) 

1500 

30.03 

117.14 

0.299* 

1600 

30.65 

119.09 

0.305* 

1700 

31.40 

121.03 

0.311 

1800 

32.28 

122.98 

0.318 

1900 

33.30 

124.92 

0.324 

2000 

34.45 

126.87 

0.330 

2100 

35.73 

128.82 

0.33  7* 

•k 

Extrapolated  values 
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TABLE  5 


Error  analysis   (at  2000  K) 


Quanti  ty 

★ 

Imprecision 

•k-k 

Inaccuracy 

Temperature 

0.5  K 

4  K 

Voltage 

0.03% 

0 . 1% 

Current 

0.03% 

0 . 1% 

Heat  capacity 

0 . 77 

370 

Electrical  resistivity 

0 . 67 

27 

Hemispherical  total 

emi  ttance 

0 . 97 

57 

Imprecision  refers  to  the  standard  deviation  of  a  quantity  as 
computed  from  the  difference  between  the  value  of  the  quantity 
and  that  from  the  smooth  function  obtained  by  the  least 
squares  method.     The  quantities  in  the  case  of  temperature, 
voltage,   and  current  are  the  individual  points  measured  in  a 
single  experiment,   and  in  the  case  of  heat  capacity,  electrical 
resistivity,   and  hemispherical  total  emittance  are  the  results 
from  all  experiments  as  tabulated  in  the  Appendix. 

** 

Inaccuracy  refers  to  the  estimated  total  error  (random  and 
systematic) . 
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3.  Discussion 

The  differences  in  the  measured  properties  for  the  three  specimens 
were  within  the  measurement  resolution  for  the  properties,   and  the  final 
smoothed  results  (represented  by  the  equations  in  Table  3  and  tabulations 
in  Table  4)  were  obtained  from  the  combined  data  for  the  three  specimens. 
The  heat  capacity,  electrical  resistivity  and  hemispherical  total 
emittance  of  zirconium  measured  in  this  work  are  presented  and  compared 
graphically  with  those  reported  in  the  literature  in  figures  I,  2  and 
3,   respectively.     The  present  results  are  for  temperatures  up  to  2100  K, 
which  is  approximately  30  K  below  the  melting  point  of  zirconium. 

The  heat  capacity  results  of  this  work  are  approximately  1%  lower 
than  those  of  Skinner  [6]  in  the  overlapping  temperature  region. 
Extrapolation  of  the  results  of  this  work  to  lower  temperatures  (1200- 
1400  K)  yields  values  which  are  3-670  lower  than  those  reported  by 
Coughlin  and  King  [7].     However,   too  much  significance  should  not  be 
attached  to  the  latter  since:     (a)  the  comparison  is  based  on  an 
extrapolation  of  100-300  K  and  (b)   the  constant  heat  capacity  reported 
by  Coughlin  and  King  [7]  is  not  realistic. 

The  electrical  resistivity  results  are  in  reasonably  good,  agreement 
(maximum  difference  less  than  1%)  with  those  of  Zhorov  [8],  and  are 
approximately  2-47Q  lower  than  those  of  Peletskii  et  al.   [9]  in  the 
overlapping  temperature  regions. 
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Zirconium  undergoes  a  solid-solid  phase  transformation  around 

1150  K.     The  measurements  of  the  geometrical  quantities  of  the  specimen 

after  a  number  of  experiments  indicated  permanent  distortions  (elongation) 

due  to  repeated  heating  and  cooling  through  the  transformation  point. 

The  reported  electrical  resistivity  results  are  corrected  for  the 

permanent  geometrical  changes .    The  magnitude  of  this  correction  was 

-8 

about  1%.    At  293  K,   the  average  electrical  resistivity  (42.8  x  10  Q-m) 

of  the  three  specimens  used  in  this  work  is  within  3.57»  of  the  values 

-8 

reported  in  the  literature  by  Adenstedt  [10]   (44.1  x  10    fi-m)  and  by  Powell 

-8  -8 

and  Tye  [11]  (42.2  x  10    Q-m  and  44.3  x  10    Q*m).     Some  of  the  differences 

in  the  electrical  resistivity  values  may  be  due  to  differences  in 
the  chemical  composition  of  the  specimens. 

The  hemispherical  total  emittance  values  reported  in  this  work 
are  higher  (10-25%)  than  those  reported  by  Timrot  and  Peletskii  [12], 
Peletskii  et  al.   [9]  and  Zhorov  [8].     Considerable  differences  in 
hemispherical  total  emittance  results  of  various  investigators  may  be 
due  to  differences  in  specimen  surface  conditions.     Changes  in  the  specimen' 
surface  conditions  were  noticed  during  this  work,  with  the  initial 
smooth  polished  surface  changing  to  an  uneven  rough  surface  as  the  result 
of  repeated  heating  and  cooling  through  the  transformation  point.  This 
may  partially  account  for  the  high  emittance  values. 
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TABLE  A -3 


Experimental  results  on  hemispherical  total  emittance  of  zirconium 


specimen-1 

i —          ■■  ■                '  " 

specimen-  2 

Temperature 

r  / 

* 

Temperature 
v  *w 

e 

* 

Ae 
(%) 

1674.7 

0.307 

1  

1 

0.87 

1650.5 

0.301 

-  2.37 

1680.1 

0.309 

0.33 

1655.5 

0.304 

1685.7 

0.310 

0.12 

1660.7 

0.307 

-  0.58 

1691.4 

0.312 

+ 

0.40 

1666.0 

0.310 

+  0.28 

1697.1 

0.313 

+ 

0.61 

1671.4 

0.312 

+  0.81 

1703.0 

0.314 

+ 

0.80 

1676.9 

0.315 

+  1.65 

1848.8 

0.320 

m 

0.24 

1980.7 

0.326 

-  0.96 

1856.9 

0.322 

+ 

0.23 

1991.0 

0.327 

-  0.86 

1865.3 

0.324 

+ 

0.68 

2001.6 

1 

0.328 

-  0.75 

1873.7 

0.325 

+ 

0.82 

2012.3 

0.329 

-  0.65 

1882.4 

0,327 

+ 

1.26 

2023.3 

0.329 

-  0.87 

1891.3 

0.328 

; 

+ 

1.39 

2034.5 

0.330 

-  0.78 

2003.3 

0.330 

0.18 

. "  i 

2014.3 

0.332 

+ 

0.22 

i 

2025.6 

0.333 

+ 

0.30 

! 
1 

I 

2037.1 

0.334 

+ 

0.38 

2048.9 

0.335 

+ 

0.45 

2061.0 

0.336  + 

0.52 

The  quantity  Ae  is  percentage  deviation  of  the  individual  results  from 
the  smooth  function  represented  by  the  pertinent  equation  in  table  3. 
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Chapter  4 


MEASUREMENT  OF  MELTING  POINT,   RADIANCE  TEMPERATURE 
(AT  MELTING  POINT),  AND  ELECTRICAL  RESISTIVITY 
(ABOVE  2100  K)  OF  ZIRCONIUM  BY  A  PULSE 

I 

HEATING  METHOD* 

A.  Cezairliyan  and  F.  Righini 
Institute  for  Materials  Research 
National  Bureau  of  Standards 
Washington,  D.  C.  20234 

Abst  ract 

A  subsecond  duration  pulse  heating  method  is  used  to  measure  the 

melting  point,   radiance  temperature  (at  650  nm)  at  the  melting  point, 

and  electrical  resistivity  (above  2100  K)  of  zirconium.     The  results 

yield  a  value  of  2128  K  for  the  melting  point  on  the  International 

Practical  Temperature  Scale  of  1968.     The  radiance  temperature  (at 

650  nm)  of  zirconium  at  its  melting  point  is  1940  K,  and  the  corresponding 

normal  spectral  emittance  is  0.367.    At  2100  K  electrical  resistivity 
-8 

is  128.7  x  10      Cl-m.     Estimated  inaccuracy  is:  8  K  in  the  melting  point 
and  in  the  radiance  temperature ,^3%  in  the  normal  spectral  emittance  and 
in  the  electrical  resistivity. 
7" 
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1 .     Int  roduc  ti on 

A  millisecond  resolution  pulse  heating  technique  was  developed 
earlier  in  connection  with  the  measurement  of  selected  thermophysical 
properties  of  electrical  conductors  at  high  temperatures  [1].  This 
technique  was  also  applied  to  the  measurement  of  the  melting  point  of 
molybdenum  [2],   tungsten  [3]  and  niobium  [41,  and  the  radiance  temp- 
erature (at  650  nm)  of  niobium  at  its  melting  point  [5]. 

In  the  present  study,   the  same  technique  is  used  for  the  measure- 
ment of  the  melting  point,   radiance  temperature  at  the  melting  point, 
and  electrical  resistivity  (above  2100  K)  of  zirconium. 

The  method  is  based  on  rapid  resistive  self-heating  of  the 
specimen  from  room  temperature  to  high  temperatures  (above  1500  K)  in 
less  than  one  second  by  the  passage  of  an  electrical  current  pulse 
through  it;  and  on  measuring,  with  millisecond  resolution,  such 
experimental  quantities  as  current  through  the  specimen,  potential  drop 
across  the  specimen,  and  specimen  temperature.     Temperature  is  measured 
with  a  high-speed  photoelectric  pyrometer  [6],  which  permits  1200 
evaluations  of  specimen  temperature  per  second.     Details  regarding  the 
construction  and  operation  of  the  measurement  system,   the  methods  of 
measuring  experimental  quantities,  and  other  pertinent  information,  such 
as  the  formulation  of  relations  for  properties,   error  analysis,  etc. 
are  given  in  earlier  publications [ 1,  7]. 
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2.  Measurements 

Measurements  were  performed  on  zirconium  specimens  in  the  form 
of  tubes  and  strips,  all  of  the  same  grade  and  obtained  from  the  same 
manufacturer.  The  specimens  were  99.9$7o  pure  and  manufacturer's  typical 
analysis  indicated  the  presence  of  the  following  impurities  in  ppm  by 
weight:  0,  125;  Hf,  40;  Fe,  30;  C,  6;  H,  3.3;  Al,  3;  N,  2.1;  Ni,  Si, 
1.5  each;  and  Ti,  1.  The  total  amount  of  all  other  detected  elements 
was  less  than  6  ppm,  each  element  being  below  1  ppm  limit. 

Melting  point  and  electrical  resistivity  measurements  were  performed 
on  three  zirconium  specimens  in  the  form  of  tubes,  designated  as 
specimen-I,  II  and  III.     The  nominal  dimensions  of  the  tubes  were: 
length,   76.2  mm;  outside  diameter,  6.3  mm  and  thickness,  0.25  mm.  A 
small  rectangular  hole  (0.5  x  1  mm)  fabricated  in  the  wall  at  the 
middle  of  the  specimen  approximated  blackbody  conditions  for  temperature 
measurements . 

The  radiance  temperature  measurements  at  the  melting  point  were 
performed  on  thirteen  zirconium  specimens  in  the  form  of  strips.  The 
nominal  dimensions  of  the  strips  were:     length,  50.8  mm;  width,  A . 7  mm 
and  thickness,   0.24  mm.    The  thin  oxide  coating  present  on  the  specimens 
(strip)  in  the  "as  received"  condition  was  removed  using  AbrdS iV2  .  Four 
different  grades  of  Abrasive,    were  used  yielding  four  different  surface 
roughnesses  (ranging  from  0.23  Mm  to  0.95  \Jxa)  for  different  specimens. 
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All  the  experiments,  with  one  exception,  were  performed  with  the 

-3  -2 

specimen  in  a  vacuum  environment  of  approximately  1.3  x  10  N-m 

(~  10      torr),  and  in  one  case  the  specimen  was  in  an  argon  environment 

at  atmospheric  pressure.     The  heating  rate  for  the  tubular  specimens 

was  approximately  3000  K-s  ^  which  corresponds  to  a  specimen  heating 

period  (from  room  temperature  to  melting  point)  of  approximately  600 

ms.    The  heating  rate  for  the  strips  varied  from  400  K-s  *  to  3200  K«s  , 

and  the  corresponding  heating  periods  changed  from  1200  ms  to  280  ms . 

The  magnitude  of  the  current  pulses  near  the  melting  point  was  about 

750  A  in  the  case  of  tubular  specimens  and  ranged  from  40A  to  300A 

for  the  9trips. 

The  high-speed  pyrometer  and  the  neutral  density  filters  used 
for  its  operation  at  different  ranges  were  calibrated  before  and  after 

the  entire  set  of  experiments,  using  a  tungsten  filament  standard  lamp, 

a  ii 

which  in  turn  was  calibrated  against  the  NBS  Temperature  Standard. 

All  temperatures  reported  in  this  paper,  except  where  explicit^, 
noted,  are  based  on  the  International  Practical  Temperature  Scale  of 
1968  [8]. 

3.     Experimental  results 
3.1.    Melting  point 

Temperature  of  the  tubular  specimens  was  measured  near  and  during 
the  initial  melting  period,  until  the  specimen  collapsed.  Typical 
results  for  the  variation  of  the  specimen  temperature  as  a  function  of 
time  (for  specimen  I)  are  shown  in  figure  1.     The  plateau  in  temperature 
indicates  the  region  of  solid  and  liquid  equilibrium.     The  complete  results 
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Table  1 

Experimental  results  of  the  melting  point  of  tubular  zirconium  specimens 


Specimen 
number 

Number  of  tempera- 
tures at  plateau 

Melting  point 

/rr\ 
(K) 

Standard  deviation 

I 

73 

2127.3 

0.4 

II 

20 

2127.5 

0.4 

III 

28 

2127.8 

0.6 
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for  the  three  tubular  airconium  specimens  are  presented  in  table  1. 

The  value  of  the  melting  point  was  obtained  by  averaging  temperature 

points  on  the  plateau  of  a  given  specimen.     The  duration  of  the 

plateau  for  specimens  II  and  III  was  shorter  than  that  of  specimen  I. 

This  may  be  attributed  to  an  early  collapse  of  these  specimens.  To 

determine  the  trend  of  measured  temperatures  at  the  plateau,  temperatures 

for  specimen  I  were  fitted  to  a  linear  function  in  time  using  the  least 

squares  method.     The  slope  of  the  linear  function  was  4.7  K-s"1,  which 

corresponds  to  a  maximum  temperature  difference  of  less  than  0.3  K  between 

the  begining  and  the  end  of  the  plateau.    This  procedure  gave  a  standard 

deviation  of  0.4  K,   the  same  as  that  obtained  by  averaging  the  temperatures. 

Values  of  the  same  order  of  magnitude  were  found  for  specimens  II  and  III. 

The  average  melting  point  of  the  three  specimens  is  2127.5  K  with  an 

from  tto.  rne^-n 

average  aksofatc.    deviation  of  0.3  K.     It  may  be  concluded  that  the 

A 

melting  point  of  zirconium  measured  in  this  work  is  2128  K. 

3.2.     Radiance  temperature  at  the  melting  point 

Radiance  temperature  measurements  were  performed  on  strips  at  650  nm 
which  corresponds  to  the  effective  wavelength  of  the  pyrometer's  inter- 
ference filter.     The  bandwidth  of  the  filter  was  10  nm.    The  circular 
area  viewed  by  the  pyrometer  was  0.2  mm  in  diameter. 

Radiance  temperature  of  zirconium  at  its  melting  point  for  the 
thirteen  experiments  (corresponding  to  thirteen  specimens)  and  other 
pertinent  results  are  reported  in  table  2.    The  variation  of  radiance 
temperature  as  a  function  of  time  near  and  at  the  melting  point  is 

shown  in  figure  2  for  three  typical  experiments  representing  three 
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specimens  with  different  initial  surface  conditions.     The  magnitude  of 
the  spike  before  the  melting  plateau  for  two  of  the  specimens  is  related 
to  the  degree  of  initial  surface  roughness  of  the  specimens.  However, 
regardless  of  the  initial  surface  and  operational  conditions,  radiance 
temperature  at  the  melting  plateau  is  approximately  the  same  for  all  the 
specimens . 

A  single  value  for  the  radiance  temperature  at  the  plateau  for  each 
specimen  was  obtained  by  averaging  the  temperatures  at  the  plateau.  The 
number  of  temperatures  used  for  averaging  ranged  from  12  to  102,  depending 
both  on  the  heating  rate  and  on  the  behavior  of  the  specimen  during  melting. 
The  standard  deviation  of  an  individual  temperature  from  the  average  was 
in  the  range  0.2  to  0.3  K  for  all  the  experiments.     Similar  values  (for 
standard  deviation)  were  obtained  when  fitting  the  temperature  data 
corresponding  to  the  premelting  period  to  a  quadratic  function  in  time. 
This  indicates  that  during  melting  no  undesirable  effects  took  place,  such 
as  vibration  of  the  specimen,   development  of  hot  spots  in  the  specimen 
and  random  changes  in  the  specimen  surface  conditions. 

To  determine  the  trend  of  measured  temperatures  at  the  plateau, 

temperatures  for  each  experiment  were  fitted  to  a  linear  function  in  time 

using  the  least  squares  method.     The  detailed  results  are  reported  in 

table  2.     The  temperature  difference  between  the  beginning  and  the  end  of 

La 

the  plateau  (corresponding  to  the  slope  in  the  plateau)  is^the  range 
0.1-0.5  K.    The  standard  deviation  of  an  individual  temperature  from 
the  linear  function  was^the  same  as  the  standard  deviation  obtained 
by  direct  averaging  of  the  temperatures. 
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The  average  radiance  temperature  at  the  melting  point  for  the 
thirteen  zirconium  specimens  was  1940.5  K  with  an  average  absolute 
difference  of  0.3  K  and  a  maximum  difference  of  0.7  K.     The  results  are 
presented  in  figure  3.     It  may  be  concluded  that  the  radiance  temperature 
(at  650  nm)  of  zirconium  at  its  melting  point  is  1940  K. 

The  normal  spectral  emittance  at  the  melting  point  was  determined 
using  the  results  of  the  radiance  temperature  (obtained  from  the  measure- 
ments on  strip-shape  specimens)  and  the  melting  point  (obtained  from  the 
measurements  on  tube-shape  specimens).     The  results  yield  a  value  of 
0.367  for  the  normal  spectral  emittance  (at  650  nm)  at  the  melting  point 
of  zirconium. 

Several  additional  experiments  were  conducted  on  zirconium  strips  in 
the  "as  received"  condition.     The  results  of  a  typical  experiment  are 
shown  in  figure  4.     The  temperature  plateau  during  melting  has  a 
considerable  slope  and  the  temperature  values  are  2  to  5  K  higher  than 
the  values  obtained  on  specimens  with  clean  surfaces.     This  may  be 
attributed  to  the  presence  of  a  thin  oxide  coating  on  the  specimens;  the 
high  temperature  values  are  consistent  with  the  general  behavior  of  the 
zirconium-oxygen  system  [9], 

3.3     Electrical  resistivity 

Electrical  resistivity  of  the  tubular  zirconium  specimens  was 

computed  using  the  relation  p  =  RA/L  where  R  is  the  resistance,  A  the 

cross-section  area  and  L  the  length  of  the  specimen  between  the  potential 

probes.     The  dimensions  were  based  on  their  room  temperature  values; 

the  cross-section  area  was  determined  from  the  measurement  of  weight  and 
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density.     Resistivity  values  above  2100  K  are  available  only  for  specimens 
II  and  III,   since  the  melting  point  determination  of  specimen  I  was  made 
without  the  voltage  probes. 

Zirconium  undergoes  a  solid-solid  phase  transformation  around 
1150  K  [10].    As  mentioned  in  an  earlier  publication  [11],   the  measure- 
ments of  the  geometrical  quantities  of  a  zirconium  specimen  after  a  number 
of  experiments  indicated  permanent  distortions  (elongation)  due  to 
repeated  heating  and  cooling  through  the  transformation  point.  The 
reported  electrical  resistivity  results  were  corrected  for  the  permanent 
geometrical  changes.    The  results  for  the  electrical  resistivity  of 
specimen  II  above  2100  K  are  shown  in  figure  5.     The  average  absolute 
difference  in  the  measured  electrical  resistivity  of  the  two  zirconium 
specimens  in  the  range  2100-2125  K  was  approximately  0.2%.     The  final 
results  for  the  electrical  resistivity  of  zirconium  were  obtained  by 
averaging  the  results  for  the  two  specimens.     This  yielded  the  values 
128.7  x  10"8  0-m  at  2100  K  and  128.9  x  10"8  fi-m  at  2120  K.     As  shown 
in  figure  5,  the  electrical  resistivity  continued  to  increase  during 
melting  of  the  specimen. 

4.     Estimate  of  errors 

Sources  and  estimates  of  errors  in  experiments  similar  to  the  ones 
conducted  in  this  study  are  given  in  detail  in  earlier  publications [ 1, 
2].     Specific  items  in  the  error  analysis  were  recomputed  whenever 
the  present  conditions  differed  from  those  in  the  earlier  publications. 
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A  summary  of  the  results  on  imprecision    and  inaccuracy    of  some  of  the 
measured  and  computed  quantities  is  given  in  table  3. 

It  may  be  seen  from  table  3  that  the  imprecision  of  blackbody 
temperature  measurements  during  heating  of  the  specimen  (before 
reaching  the  melting  point)  was  approximately  the  same  as  that  during 
the  melting  period.     This  indicates  that  in  the  experiment  the  initial 
melting  phase  progressed  normally  and  that  there  were  no  undesirable 
effects  during  melting  of  the  specimen,   such  as  vibration  of  the  specimen, 
movement  of  the  blackbody  sighting  hole  in  the  specimen,  and  the 
instantaneous  development  of  hot  spots  or  zones  in  the  specimen.  Also 
the  imprecision  of  the  radiance  temperature  measurements  was  the  same 
as  that  of  blackbody  temperature  measurements,  which  indicates  that 
the  surface  conditions  of  the  specimen  did  not  change  in  a  random  fashion 
during  the  melting  period. 

5.  Discussion 

The  values  of  the  melting  point  of  zirconium  reported  in  the 
literature  are  given  in  table  4.     No  large  discrepancies  are  found:  all 
values  are  within  the  combined  inaccuracies  of  the  measurements.  There 

^"Imprecision  refers  to  the  standard  deviation  of  an  individual  point  as 
computed  from  the  difference  between  measured  value  and  that  from  the 
smooth  function  obtained  by  the  least  squares  method. 

2 

Inaccuracy  refers  to  the  estimated  total  error  (random  and  systematic) . 
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Table  3 


Imprecision  and  inaccuracy  of  measured  and  computed  quantities 


Quantity 

Imprecision 

Inaccuracy 

Temperature:  blackbody 

0.5  K 

5  K 

(during  heating) 

Temperature:     blackbody  4 

0.5  K 

5  K 

(during  melting) 

Temperature:  radiance 

0.3  K 

5  K 

(during  melting) 

Voltage 

V  .  U  J  u 

U  .  i/o 

Current 

0.03% 

0.1% 

Melting  point  of  zirconium 

0.5  K 

8  K 

Radiance  temperature  of  zirconium 

0.3  K 

8  K 

at  its  melting  point 

Normal  spectral  emittance 

0.3% 

3%  I 

Electrical  resistivity 

0.05% 

3% 
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are  indications  that  the  specimens  used  in  some  of  the  earlier  measure- 
ments were  less  pure  than  those  used  in  the  present  study.     Some  of  the 
high  values  may  be  attributed  to  relatively  large  hafnium  (melting 
point  about  2500  K)  content  of  the  specimens. 

Normal  spectral  emittance  (at  0.645  nm)  of  zirconium  at  its  melting 
point  is  reported  to  be  0.318  [ 19 ] ;  this  is  approximately  13%  lower 
than  the  value  obtained  in  this  work.     Since  in  both  cases,  the  same 
value  for  the  melting  point  was  considered,   the  difference  in  emittance 
may  be  attributed  to  the  difference  in  radiance  temperature  at  the 
melting  point,  which  is  22  K.     No  satisfactory  explanation  for  this 
large  discrepancy  has  been  found. 

In  an  earlier  publication  [11],  measurements  of  the  electrical 

resistivity  of  zirconium  were  reported  for  temperatures  up  to  2100  K.  At 

2100  K  the  difference  between  the  present  value  and  that  reported  earlier 

is  approximately  0.17o.     Electrical  resistivity  of  zirconium  up  to  2000  K 

were  reported  by  Peletskii  et  al.   [20]  and  by  Zhorov  [21].  Extrapolation 

-8 

of  those  results  to  2100  K  give  the  values  127.9  x  10      Q-m  and 
-8 

132.0  x  10      Q-m,   respectively.     The  results  of  this  work  differ  from 
the  extrapolated  values  by  0.7%  in  the  case  of  Peletskii  et  al.  [20] 
and  by  2 . 6%  in  the  case  of  Zhorov  [21].  Electrical  resistivity  behaved 
normally  until  approximately  5-10  K  below  the  melting  point  (figure  5). 
Above  this,  a  departure  from    normalcy    was  observed  which  may  be  due 
to:     (a)  premelting  effects  of  impurities  present  in  the  specimen,  and 
(b)     increase  in  vacancy  concentration.     Small  temperature  gradients 
in  the  specimen  also  may  partially  account  for  this. 
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As  it  may  be  noted  from  table  2,   the  linear  fitting  of  radiance 
temperatures  at  the  plateau  shows  a  bias  toward  positive  slopes.  Two 
possible  explanations  may  be  given  for  this  gradual  increase:  (a) 
a  small  change  in  the  normal  spectral  emittance  of  zirconium  during 
melting,  and  (b)  partial  formation  of  zirconium  oxide  on  the  specimen 
surface,  due  to  the  high  chemical  reactivity  of  zirconium.     It  may  be 
noted  that,  with  the  present  system,  it  was  not  possible  to  follow  the 
entire  melting  process  because  the  specimen  collapsed  and  opened  the  main 
electrical  circuit  prior  to  the  completion  of  melting.  Approximate 
calculations  in  the  case  of  zirconium  strips  indicate  that  approximately 
one  third  of  the  specimen  was  melted  during  the  experiment. 

In  conclusion,  the  results  of  this  study  have  yielded  new  values 
for  the  melting  point  and  the  electrical  resistivity  of  zirconium  above 
2100  K,  and  have  shown  the  constancy  and  reproducibility  of  the 
radiance  temperature  of  zirconium  at  its  melting  point.     The  scheme 
of  measuring  radiance  temperature  at  the  melting  point,  so  far 
demonstrated  for  niobium  [5]  and  zirconium,  may  be  an  easy,  accurate 
and  practical  way  for  performing  secondary  calibrations  and  for 
conducting  overall  on-the-spot  checks  on  complicated  measurements 
systems  at  high  temperatures. 
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Chapter  5 


Thermodynamic  Studies  of  the 
or  -»  6  Phase  Transformation  in  Zirconium  Using 

•k 

a  Subsecond  Pulse  Heating  Technique 


•irk 

A.  Cezairliyan  and  F.  Righini 
Institute  for  Materials  Research 
National  Bureau  of  Standards 
Washington,  D.  C.  20234 

Abstract 

Measurements  of  the  temperature  and  energy  of  the  oe  -*  3  phase 
transformation,  and  the  electrical  resistivity  near  and  at  the 
transformation  point  of  zirconium  using  a  subsecond  duration  pulse 
heating  technique  are  described.     The  results  yield  1147  K  for  the 
transformation  temperature  and  3980  J-mol  *  for  the  transformation 
energy.     Electrical  resistivity  is  found  to  decrease  by  177o  during  the 
transformation.     Estimated  inaccuracies  of  the  measured  properties  are: 
10  K  for  ^transformation  temperature,  57Q  for  ^transformation  energy,  and 
2%  f or^electrical  resistivity. 
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1.  Introduction 

In    3.no.ti\t2<c  publication  [1]*,   the  applicability  of  a  rapid 
pulse  heating  technique  to  studies  of  solid-solid  phase  transformations 
at  high  temperatures  was  demonstrated  by  measurements  on  iron  at  the 
X  ~*  6   transformation  point. 

The  objective  of  this  paper  is  to  apply  the  same  technique  to 
studies  of  the  a  -»  p  phase  transformation  in  zirconium.    Measurements  of 

the. 

the  temperature  and  energy  ofphase  transformation,  and  electrical 
resistivity  near  and  at  the  transformation  point  are  reported. 

The  method  is  based  on  rapid  resistive  self-heating  of  the  specimen 
from  room  temperature  to  high  temperatures  in  less  than  one  second  by 
the  passage  of  an  electrical  current  pulse  through  it;  and  on  measuring, 
such  experimental  quantities  as  current  through  the  specimen,  potential 
drop  across  the  specimen,  and  specimen  temperature.    Recordings  of  the 
experimental  quantities  are  made  digitally  every  0.4  ms  with  a  full- 
scale  signal  resolution  of  approximately  one  part  in  8000.  Specimen 
temperature  is  measured  with  a  high-speed  photoelectric  pyrometer  [2]. 
The  details  regarding  the  construction  and  operation  of  the  measurement 
system  are  given  in  earlier  publications  [3,  4], 


Figures  in  brackets  indicate  the  literature  references  at  the  end 
of  this  paper. 
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2.  Measurements 

The  measurements  were  made  on  three  zirconium  specimens  of  99,98% 
purity.    The  specimens  were  tubes  fabricated  from  rods  by  removing  the 
center  portion  using  an  electro-erosion  technique.    The  nominal  dimensions 
of  the  specimens  were:     length,  76.2  mm;  outside  diameter,  6.3  mm;  and 
wall  thickness,  0.5  mm.    The  outer  surfaces  of  the  specimens  were 
polished  to  reduce  heat  loss  due  to  thermal  radiation.    According  to 
the  manufacturer's  analysis,  the  specimens  contained  the  following 
impurities,  in  ppm  by  weight:     0,   125;  Hf,  40;  Fe,  30;  C,  6;  H,  3.3; 
Al,  3;  N,  2.1;  Ni  and  Si,  1.5  each;  and  Ti,  1.    The  total  amount  of 
all  other  detected  elements  was  less  than  6  ppm,  each  element  being 
below  1  ppm  limit. 

Duration  of  the  current  pulses  ranged  from  220  to  620  ms.  Specimen 
heating  rates  varied  from  1800  Ks  *  to  4000  Ks         All  the  experiments 
were  conducted  with  the  specimens  in  a  vacuum  environment  of  approximately 
1.3  x  10"3  N-nf2  (~  10"5  torr). 

The  a  -»  p  transformation  (from  hexagonal  close-packed  form  to  body- 
centered  cubic)  was  manifested  by  a  plateau  in  the  temperature  versus 
time  relation  for  the  specimen  during  heating.    The  transformation 
temperature  for  each  specimen  was  obtained  by  averaging  the  temperatures 
at  the  plateau. 
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The  transformation  energy  was  obtained  from  the  time  integral  of 
the  power  absorbed  by  the  specimen  during  the  transformation  as 
defined  by  the  plateau.    The  instantaneous  value  of  the  absorbed  power 
was  obtained  by    subtracting  the  power  loss  due  to  thermal  radiation  from 
the  imparted  power  (current  through  the  specimen  times  the  potential 
drop  across  the  specimen).    The  hemispherical  total  emittance  needed 
at  1147  K  was  obtained  from  the  extrapolation  of  the  measurements 
above  1500  K  reported  in  an  earlier  publication  [5].    Since  at  the 
transformation  temperature  power  loss  from  the  specimen  due  to  thermal 
radiation  was  approximately  0.7%  of  the  imparted  power,  even  considerable 
uncertainties  (10%)  in  the  emittance  do  not  contribute  any  significant 
errors  to  the  absorbed  power. 

3.     Experimental  results 

Except  where  explicitly  noted,  all  temperatures  reported  in  this 
paper  are  based  on  the  International  Practical  Temperature  Scale  of 
1968  [6]r    In  all  computations,  the  geometrical  quantities  are  based  on 
their  dimensions  at  room  temperature  (298  K) .     The  experimental  results 
for  the  temperature  and  energy  of  the  a  -».f3  phase  transformation  are 
presented  in  table  1. 

The  average  value  for  the  transformation  temperature  is  1146.6  K 
with  a  maximum  and  average  absolute  difference  from  this  value  of    ' 8  K 
and  1.1  K,  respectively.    The  average  value  for  the  transformation 
energy  is  3977  J-mol  ^  with  a  maximum  and  average  absolute  difference 
from  this  value  of  0.4%  and  0.3%,  respectively.     It  may  be  concluded 

that  for  ot  -*  (3  transformation  in  zirconium,  transformation  temperature 
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is  1147  K  and  the  transformation  energy  is  3980  J*mol 

As  discussed  in  an  earlier  publication  [5],  the  measurements  of 
the  geometrical  quantities  of  a  zirconium  specimen  after  a  number  of 
experiments  indicated  permanent  distortions  (elongation)  due  to 
repeated  heating  and  cooling  through  the  transformation  point.  The 
values  reported  in  this  paper  correspond  to  the  first  experiment  for 
a  specimen  for  which  no  correction  was  needed. 

Typical  experimental  results  are  presented  in  figures  1,  2,  and 
3,  all  referring  to  specimen  3.  Specimen 

temperature  near  and  at  the  transformation  point    is    shown  in  figure  ). 

Figure  2  shows  the  variation  of  electrical 

resistivity  as  a  function  of  time,  while  the  variation  of  resistivity 

with  temperature  is  shown  in  figure  3.    The  dashed  lines  in  figure  3  are 

the  best  least  squares  fittings  of  the  resistivity  versus  temperature 

data  in  the  respective  regions.     It  may  be  seen  that  the  change  in  the 

electrical  resistivity  (100(0,    _      -  p  -      )  during  the  transformation 

■  before  after 

was  about  1770.  Extrapolation  of  the  resistivity  values  about  the 
transformation  to  1500  K  is  in  agreement,  within  0.37o,  with  those 
reported  in  an  earlier  publication  [5]. 

The  details  for  estimating  errors  in  measured  and  computed  quantities 
using  the  present  measurement  system  are  given  in  an  earlier  publication 
[4].    In  this  paper  the  specific  items  were  recomputed  whenever  the 
present  conditions  differed  from  those  in  the  earlier  publication. 
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The  results  for  imprecision  and  inaccuracy  are:  3  K  and  10  K 
for  the  transformation  temperature,  17,  and  5%  for  the  transformation 
energy,  and  0.57,  and  27,  for  the  electrical  resistivity. 

In  the  case  of  temperature,    the  standard  deviation  of  an 
individual  point  from  the  average  transformation  temperature  in  a  given 
experiment  is  approximately  3  K.     For  the  three  experiments,  the 
average  absolute  difference  of  transformation  temperatures  from  their 
mean  is  1.2  K.    In  the  case  of  energy,   the  imprecision  includes  the 
uncertainty  in  the  computation  of  the  time  duration  of  the  transformation 
plateau.     Different  ways  in  selecting  temperature  points  belonging  to 
heating  or  plateau  regions  lead  to  average  energy  differences  of  0.67o. 

4.  Discussion 

The  first  evidence  of  a  possible  solid-solid  phase  transformation 
in  zirconium  was  presented  by  Zwikker  in  1926  [71 .     Since  his  experiments 
were  performed  on  zirconium  heated  in  air,   the  transformation  occurred 
over  an  extended  temperature  range.    A  sharp  transformation  was  reported 
by  Vogel  and  Tonn  [81  from  dilatometric  and  metallographic  measurements 
and  was  confirmed  by  de  Boer  et  al.   [91,  who  noted  that  the  temperature 

2 

Imprecision  refers  to  the  standard  deviation  of  an  individual  point  as 
computed  from  the  difference  between  measured  value  and  either  the 
average  or  the  smooth  function  obtained  by  the  least  squares  method. 

3 

Inaccuracy  refers  to  the  estimated  total  error  (random  and  systematic). 
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range  was  extended  by  small  amounts  of  oxygen  and  nitrogen  in  zirconium. 
The  values  for  the  a  -*  Q  transformation  temperature  for  zirconium 
reported  in  the  literature  are  given  in  table  2.     Some  of  the  old 
results  could  not  be  corrected  to  IPTS-68,  due  to  insufficient 
information  in  the  original  papers.    However,   the  contribution  of  the 
possible  differences  in  the  temperature  scales  may  not  exceed  a  few 
degrees  at  the  transformation  temperature. 

The  present  value  for  the  transformation  temperature  is  higher 
than  most  of  the  values  given  in  the  literature.     This  may  partly  be 
due  to  the  high  heating  rates  used  in  this  work.    A  similar  effect, 
though  in  the  opposite  direction,  was  observed  by  Duwez  [10]  in 
experiments  on  rapid  quenching  of  zirconium  specimens  (cooling  rates 
from  4  to  10,000  K«s  The  transformation  temperature  was  found 

to  decrease  with  increasing  cooling  rate  by  as  much  as  15  K.  This 
behavior  was  confirmed  by  Hayes  and  Kaufman  [15].     The    result  of  the 
energy  of  transformation  of  zirconium  obtained  in  this  work  (3980  J«mol~*) 
is  67,  higher  than  the  value  (3740  J.mol"1)  of  Douglas  and  Victor  [14], 
37c  higher  than  the  value  (3850  J.mol"1)  of  Coughlin  and  King  [16],  and 
0.17„  higher  than  the  most  recent  value  (3975  J.mol"1)  of  Vollmer  et  al. 
[13]. 

The  general  trend  of  the  behavior  of  electrical  resistivity  near 
and  at  the  transformation  point  of  zirconium  is  as  follows  (figures  2. 
and  3).     Resistivity  reaches  a  maximum  about  10  K  below  the  transformation 
point,   then  it  decreases  slightly  and  undergoes  a  large  change  (decrease) 
at  the  transformation  point,  after  which  it  decreases  some  more  and 
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Table  2 


s  of  the  ot-Q  transformation  temperature  of  zirconium  reported  in  the  literatu 


Investigator 

Reference 

Year 

Temperature3 
(K) 

Vogel  and  Tonn 

o 
o 

1  Q1  1 

iyj  l 

1135  ±  5 

de  Boer  et  al . 

Q 

y 

LyjO 

1  no    _i_    i  r\ 
LLjo  ±  1U 

Duwez 

IU 

1  QC1 

iy  j  l 

ll^J  to  llJo 

Domagala  and  McPherson 

1 1 
1 1 

lyD't 

1  IOC 

Kneip  and  Betterton 

12 

1956 

1138  to  1146 

Vollmer  et  al. 

13 

1967 

1155 

Present  work 

1147  ±  10 

All  temperatures  in  this  table  are  those  reported  by  the  authors  of  the 
original  papers.     No  correction  to  IPTS-68  was  possible,  due  to 
insufficient  information  in  the  original  papers. 
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reaches  a  minimum  at  about  5-10  K  above  the  transformation  point. 
Above  this  minimum,  resistivity  increases  with  temperature  almost 
linearly . 

Because  of  the  strong  dependence  of  the  electrical  resistivity 
on  composition  for  zirconium  at  the  transformation  point,   it  is 
difficult  to  compare  resistivity  values  reported  in  the  literature  on 
an  absolute  basis.     However,  a  meaningful  comparison  may  be  made  based 
on  the  ratio  of  the  maximum  to  the  minimum  resistivities  for  a  given 
specimen.     The  results  of  such  a  comparison  are  given  in  table  3. 
High  values  for  the  ratio  and  small  temperature  range  during  which 
transformation  is  completed  are  indications  of  the  high  purity  of 
the  specimens. 
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Chapter  6 


RADIANCE  TEMPERATURE  (at  650  nm)  OF  MOLYBDENUM 
AT  ITS  MELTING  POINT 

Ared  Cezairliyan 
National  Bureau  of  Standards 
Washington,  D.  C.  20234 


Measurements  of  the  radiance  temperature  (at  650  nm)  of  molybdenum 
(99.95  percent  pure)  at  its  melting  point  were  performed  on  14 
specimens  in  the  form  of  strips.     The  nominal  dimensions  of  the  strips 
were:     length,   50.8  mm;  width,  4 . 7  mm  and  thickness,   0.24  mm.  Before 
the  experiments,   specimen  surface  was  treated  using  abrasive.  ■  three 
different  grades  of  c2br<3Sii/£    were  used  yielding  three  different 
surface  roughnesses   (ranging  from  approximately  0.2  to  0.5  Mm)  for 
different  specimens. 

All  the  experiments  were  performed  with  the  specimen  in  an  argon 
environment  at  atmospheric  pressure.     The  heating   rates  for  different 
specimens  were  in  the  range  1200  K«s  *  to  3500  K*s     ,  corresponding  to 
specimen  heating  periods  (from  room  temperature  to  the  melting  point) 
in  the  range  0.4  -  1  s. 
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Radiance  temperature  of  molybdenum  at  its  melting  point  for  the 
14  specimens  and  other  pertinent  results  are  reported  in  Table  1. 
The  variation  of  radiance  temperature  as  a  function  of  time  near  and 
at  the  melting  point  is  shown  in  Figure  1  for  three  typical  experiments 
representing  three  specimens  with  different  initial  surface  conditions. 
The  magnitude  of  the  spike  before  the  melting  plateau  for  the  specimens 
is  related  to  the  degree  of  initial  surface  roughness  of  the  specimens. 
However,   regardless  of  the  initial  surface  and  operational  conditions, 
radiance  temperature  at  the  melting  plateau  is  approximately  the  same 
for  all  the  specimens.     Temperature  of  the  spike  in  all  the  experiments 
ranged  from  4  to  16  K  above  the  plateau  temperature. 

A  single  value  for  the  radiance  temperature  at  the  plateau  for 
each  specimen  was  obtained  by  averaging  the  measured  values  of  the  temperatures 
at  the  plateau.     The  number  of  measurements  used  for  averaging  ranged  from 
25  to  99,  depending  both  on  the  heating  rate  and  on  the  behavior  of  the 
specimen  during  melting.     The  standard  deviation  of  an  individual  measurement 

from  the  average  was  in  the  range  0.3  to  0.4  K  for  all  the  experiments. 
Similar  values  (for  standard  deviation)  were  obtained  when  fitting  the 
temperature  data  corresponding  to  the  premelting  period  to  a  quadratic 
function  in  time.     This  indicates  that  during  melting  no  undesirable 
effects  took  place,   such  as  vibration  of  the  specimen,  development 
of  hot  spots  in  the  specimen  and  random  changes  in  the  specimen 
surface  conditions. 
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To  determine  the  trend  of  measured  temperatures  at  the  plateau, 
temperatures  for  each  experiment  were  fitted  to  a  linear  function  in 
time  using  the  least  squares  method.     The  detailed  results  are  reported 
in  Table  1.     The  temperature  difference  between  the  beginning  and  the 
end  of  the  plateau  (corresponding  to  the  slope  in  the  plateau)  is  in 
the  range  0  -  0.5  K.     The  standard  deviation  of  an  individual  temp- 
erature from  the  linear  function  was  approximately  the  same  as  the 
standard  deviation  obtained  by  direct  averaging  of  the  temperatures. 

The  average  radiance  temperature  at  the  melting  point  for  the 
14  molybdenum  specimens  was  2530.4  K  with  a  standard  deviation  of  0.6  K 
and  a  maximum  absolute  difference  of  1.2  K.     The  results  are  presented 
in  Figure  2.     It  may  be  noted  that  the  maximum  absolute  difference 
of  the  results  for  the  specimens  with  condition  A  surfaces  (least 
rough)  is  only  0.7  K.     The  two  experiments   (corresponding  to  specimen-2 
and  specimen-10)  with  relatively  large  negative  deviations  (about 
-1.1  K  and  -1.2  K)  have  caused  a  downward  shift  of  the  average  by 
about  0.2  K.     In  order  to  reduce  the  contribution  of  these  two 
experiments,    the  average  temperature  is  rounded-off  to  the  higher  value. 
This  is  also  compatible  with  the  average   (2530.6  K)   of  the  results  for 
the  specimens  with  condition  A  surfaces.     It  may  be  concluded  that 
the  radiance  temperature  of  molybdenum  at  its  melting  point  is  2531  K, 
with  an  estimated  inaccuracy  of  10  K. 
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The  radiance  temperature  (2531  K)  and  the  melting  point  (2894  K) 
measured  with  the  present  system  [1]  yield  a  value  of  0.334  for  the 
normal  spectral  emittance  (at  650  nm)  of  molybdenum  at  its  melting 
point . 

Several  additional  experiments  were  conducted  on  molybdenum  in 
"as  received"  condition.     The  results  of  these  experiments  were 
generally  3-5  K  higher  than  the  value  2531  K.     The  higher  values  might 
be  attributed  to  contamination  on  the  specimen  surface.     Because  of 
the  high  vapor  pressure  of  molybdenum  at  its  melting  point,  it  was 
not  possible  to  conduct  meaningful  experiments  with  the  specimen  in 
a  vacuum  environment. 

Bonnell  et  al.   [2]  have  reported  a  value  of  2510  K  for  the 
radiance  temperature  (at  645  nm)  of  molybdenum  (99.9  percent  pure) 
at  its  melting  point.     No  satisfactory  explanation  for  the  large 
discrepancy  (21  K)  between  the  result  of  Bonnell  et  al.  and  that  of 
the  present  work  has  been  found.     In  the  past  [3,  41,  differences  of 
the  same  magnitude  and  direction  were  also  noted  in  the  results  obtainpH 

bv  the  two  laboratories  for  other  refractorv  metals.  The 
differences  were  20  K  for  niobium  (at  2425  K)  and  22  K  for  zirconium 
(at  1940  K). 
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Matrix  Isolation  Studies  on  the  Reaction  of 
Iron  Vapor  with  Oxygen 

S.  Abramowitz  and  N.  Acquista 
National  Bureau  of  Standards 
Institute  for  Materials  Research 
Washington,  D.  C.  20234 

Introduction 

The  reactions  of  metallic  vapor  beams  with  oxygen  species  such  as 
N2O,  0,  O2,  and  O3,  are  often  exothermic  enough  to  populate  electronically 
excited  levels  of  the  diatomic  oxides  formed.     Recently  several  groups 
(1,2)  have  shown  that  reactions  of  Ba(g)  with  N20  and  O3  give  high  yields 
of  electronically  excited  BaO  (on  the  order  of  10%  of  the  barium  vapor 
consumed  results  in  photons  emitted  in  the  A-X  system  of  BaO) .     In  order 
to  fully  understand  such  processes  and  to  evaluate  them  for  possible  use 
as  chemical  laser  systems  a  knowledge  of  the  rates  of  all  the  reactions 
involved  and  the  nature  of  any  chemical  intermediates  is  necessary. 

The  methods  of  matrix  isolation  spectroscopy  are  particularly  useful 
for  identifying  possible  intermediates  in  oxidation  processes.  The 
infrared  spectra  of  matrix  isolated  Ba02  formed  by  the  reaction  of  Ba(g) 
with  O2  has  been  observed  (3).     Simple  oxides  of  metals  such  as  Pb(4), 
Sn(5),  Ge(6),  Ni(7),  Pt(8),  Pd(8),  Cu(9),  Cr(9),  Li(lO.ll),  Na(12),K(13) 
and  Rb(13)  have  been  identified  using  this  method.     In  addition  to  these 
species  reaction  of  alkali  metals  (14,15,16)  and  some  alkaline-earth  metals 
with  03/Ar  (17)  have  given  infrared  spectra  of  matrix  isolated  MO3  species. 
In  some  cases  where  the  AH  of  reaction    is  sufficient  to  yield  the  monoxide 
(when  the  metal  vapor  is  reacted  with  02)  the  monoxide,  dioxide  and  higher 
oxides  are  observed.     Mole  ratios  of  the  observed  oxides  can  often  be  con- 
trolled by  the  02/Ar  (or  other  matrix  gas)  used.     Such  was  the  case  for  the 
reaction  of  U(g)   (18)  and  Th(g)   (19)  with  02  or  N20  where  varying  ratios  of 
U0,  U02  and  U03  (and  ThO,  Th02)  are  observed  as  a  function  of  the  02  or  N20 
concentration.     The  U02 ,  UO3  and  Th02  formed  in  these  reactions  have  the  same 

structures  as  those  obtained  by  generating  the  metallic  oxides  using  standard 
Knudsen  effusion  techniques.     The  reaction  of  Fe(g)  with  02  and  N20  is  reported 
in  this  study.     The  reaction  of  Fe(g)  with  02  to  produce  FeO  is  endothermic  by 
about  20  kcal/mole  and  is  not  expected  to  proceed.     Rather  a  Fe02  species  is 
expected.     The  reaction  of  Fe  +  N20  to  give  FeO  +  N2  is  exothermic  by  60 
kcal/mole  and  is  expected  to  proceed. 

Experimental 

The  experimental  apparatus  used  in  these  studies  has  been  described  (18). 
An  Air  Products  closed  cycle  helium  refrigerator  was  used  to  obtain  liquid 
hydrogen  temperatures.     Temperatures  of  the  deposition  holder  were  measured 
using  a  hydrogen  bulb  thermometer  and  a  gold-iron  vs.   chromel  thermocouple. 
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A  heater  at  the  deposition  plate  allowed  temperatures  to  be  varied  to 
obtain  good  matrices  and  to  carry  out  controlled  diffusion  studies.  The 
iron  was  vaporized  from  tungsten  and  molybdenum  Knudsen  cells.     In  some 
of  the  experiments  boron  or  carbon  liners  were  used  to  minimize  reactions 
of  the  iron  with  the  tungsten  and  molybdenum.     Temperatures  of  the  Knudsen 
cell  were  measured  using  W  25%  Re  vs  W3%  Re  or  Pt  vs.  Pt/Rh  thermocouples, 
or  in  some  cases  an  optical  pyrometer.    An  electron  beam  furnace  is  used 
to  generate  the  iron  vapor  in  the  one  to  ten  micron  pressure  range.  Argon 
and  nitrogen  were  used  as  the  matrix  gases.     Isotopically  enriched  oxygen 
samples  obtained  from  the  Miles-Yeda  company  were  used  in  some  of  the 
studies . 

Infrared  spectra  were  observed  with  a  Perkin-Elmer  model  99G  mono- 
chromator  equipped  with  either  a  40  or  1502./mm  diffraction  grating  for 
the  500  and  1000  cm-1  region.     Interference  filters  were  used  to  block 
out  higher  order  radiation.     A  Perkin-Elmer  model  301  was  used  to  obtain 
infrared  spectra  in  those  parts  of  the  500  cm-1  region  where  water  vapor 
lines  are  strong.     Calibration  of  both  infrared  instruments  is  accomplished 
with  water  and  ammonia  lines.     It  is  estimated  that  the  frequencies  are 
good  to  ±0.5  cm-1  in  both  regions  studied  with  the  99G .      The  PE  301  spectra 
observed  in  the  500  cm-1  are  probably  accurate  to  +0.8  cm"1. 

Results 

The  infrared  spectrum  of  the  product  of  the  reaction  of  Fe(g)  with 
02/Ar  (1/100)  is  shown  in  Figure  1.  The  extent  of  reaction  and/or  the 
absorptivity  coefficient  of  the  species  formed  is  small  since  many  hours 
of  reaction  are  required  in  order  to  obtain  the  spectrum  shown.  Spectrum 
a  is  that  of  Fe(g)  +  1602/Ar  (1/100).  The  main  peak  in  the  spectrum  is  at 
945.9  cm-1.  Other  peaks  at  954.7,  969.5,  918.8  and  915.3  era"1  are  due  to 
polymeric  features  (or  those  possibly  involving  one  iron  atom  and  more  than 
one  02  unit).-  These  features  are  increased  in  intensity  relative  to  the 
945.9  cm-1  which  diminishes  in  intensity  when  controlled  diffusion  of 
the  matrix  at  temperatures  of  about  30K  is  allowed  to  proceed.  These  fea- 
tures are  also  more  prominent  when  lower  dilutions  of  02/Ar  are  used.  The 
945.9  cm-1  is  extremely  narrow  with  a  Avj/p  (on  an  absorbance  scale)  of 
about  1  cm  1 .  In  addition  to  this  feature  at  945.9  cm  1  another  band 
attributable  to  a  low  molecular  weight  species  appears  at  517.1  cm-1.  This 
feature  also  loses  intensity  at  about  the  same  rate  as  the  945.9  cm-1  upon 
controlled  diffusion.  Other  weak  features  in  the  300  cm-1  region  are  not 
reproducibly  present  and  do  not  behave  in  a  similar  manner  upon  controlled 
diffusion.  The  517.1  cm-1  feature  is  much  weaker  than  its  945.9  cm-1  counter- 
part. A  rough  measurement  of  peak  intensities  indicate  about  a  seven  to  one 
ratio.  The  517.1  cm"1  feature  appears  to  be  also  narrow,  however  due  to  its 
weakness  it  is  difficult  to  measure  its  half  width. 


An  isotopically  enriched  sample  of  02  was  used  in  the  spectrum  b.  The 
oxygen  contained  30.5%  1802,  19.5%  1602,  48.7%  160180  and  insignificant 
amounts  of  1 702 ,   170180,  and  160170.     The  single  peak  at  945.9  cm-1  is 
replaced  with  three  peaks  in  the  peak  intensity  ratio  of  about  2:5:3  at 
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945.9,  930.8,  and  911.4  cm-1.     (These  spectra  are  direct  reproductions 
of  those  obtained  on  the  PE  99G}     which  is  a  single  beam  instrument. 
Therefore,  the  Io  is  not  flat,  but  rather  tends  to  increase  in  going  to 
higher  frequency  due  to  the  blackbody  characteristics  of  the  Globar 
source,  transmission  of  the  interference  filter,  and  the  diffraction 
grating  efficiency.)     This  ratio  is  in  good  agreement  with  that  of  the 
isotopically  enriched  oxygen  sample  used.     In  the  500  cm-1  region  a 
feature  is  found  at  508.1  cm-* .     The  weakness  of  this  feature  precludes 
the  observation  of  the  517.1  cm-1  feature  in  the  same  spectrum.  Simi- 
larly it  wasn't  possible  to  observe  the  feature  which  will  be  seen  to 
be  assignable  to  Fe1802  in  this  spectrum. 

In  order  to  observe  the  features  due  to  Fe1802  it  was  necessary  to 
use  an  isotopically  enriched  *802  sample  con  taining  86.0%  1 802 ,  12.2% 
160180,  0.4%  1602,  and  similarly  insignificant  amounts  of  170*80,  170160 
and  an  170160    and  an  1702  concentration  about  two  orders  of  magnitude 
smaller.     The  spectrum  in  Figure  lc  is  that  of  the  product  of  the  reaction 
of  Fe(g)  with  this  oxygen  sample  diluted  in  50  parts  argon.     The  911.4  cm-1 
feature  assigned  to  Fe*802  is  prominent  in  this  spectrum.     The  930.8  cm-1 
feature,  assignable  to  Fe*60180  is  also  clearly  seen.     The  intensity  of 
this  feature  relative  to  that  of  Fe1802  is  approximately  7:1  in  agreement 
with  the  isotopic  abundances  of  1 802  and  160*80  in  the  isotopically  enriched 
oxygen  sample  used.     A  peak  at  494.0  cm-1  was  also  found.     As  was  the  case 
in  the  experiments  with  the  sample  containing  48.7%  160180  and  lesser  amounts 
of  the  other  isotopic  varieties  it  was  not  possible  to  see  any  features  due 
to  160180  in  this  region.     The  low  intensity  of  all  the  features  in  the 
500  cm-*  is  responsible  for  this.     It  might  be  added  for  completeness  that 
these  low  intensities  only  allowed  controlled  diffusion  experiments  to  be 
performed  with  the  1602  sample.     (This  is  the  only  sample  containing  essen- 
tially 100%  of  one  isotopic  variety  of  oxygen.) 

In  a  separate  set  of  experiments  the  reaction  of  Fe(g)  with  N20/Ar  has 
been  studied.  FeO  is  the  principal  product  of  this  reaction  and  was  found 
in  the  argon  matrix  at  about  880  cm-*  in  good  agreement  with  the  gas  phase 
value  derived  from  its  electronic  spectrum.  Another  peak  located  at  some- 
what higher  frequency  is  also  observed.  It  is  believed  this  peak  is  assign- 
able to  a  Fe-nitrogen  complex,  possibly  FeN2 .  However  a  proof  of  this  and 
elucidation  of  the  bonding  obviously  requires  isotopic  studies;  utilizing 
N2  and  14N15N. 

Discussion 

A  M02  species  can  either  have  a  side-on  or  end-on  structure. 

M-0=0  M 

0=0 

end-on  .  , 

side-on 

The  side-on  structure  shown  is  a  monodentate  model  while  an  isosceles 
triangle  structure  is 
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also  possible.     Both  limiting  models  of  the  side-on  structure  have 
structure.     The  end-on  structure  will  either  have  C      or  C  symmetry 
depending  upon  whether  the  molecule  is  linear  or  ben^ .     Enl-on  structures 
have  been  found  for  several  metal  nitrogen  complexes.     The  spectrum 
expected  for  the  end-on  and  both  limiting  cases  of  the  side-on  models 
in  the  0-0  stretching  region  are  shown  below  for  a  1:2:1  mixture  of 
M1602,  M160180  and  M*802. 


Mlb0ia0  end-on 
side-on 

A  three  line  spectrum  with  an  intensity  ratio  of  1:2:1  is  expected  for 
a  side-on  model.     The  central  component  shows  the  spectral  equivalence 
of  the  mixed  isotopic  species  for  this  model.     The  end-on  model  gives  a 
four  line  spectrum  with  a  closely  spaced  doublet  with  absorption  intensi- 
ties in  the  ratio  of  1:1:1:1.     The  inequivalence  but  equal  probability  of 
M180180  and  M180160  results  in  a  small  separation  and  equal  intensities 
of  the  central  doublet. 


The  MO  stretching  region  would  show  either  a  triplet  or  quartet  for 
side-on  and  end-on  models  of  MO2  respectively.     These  modes  are  expected 
to  be  as  shown  below  for  an  isotopic  abundance  of  1^02 :  2*60*  80:  11809. 


M1602 


M1  80, 


M1602  M160180 


M180160  M1802 


M^O1  °0 
side-on 


end-on 


The  side-on  model  would  exhibit  the  same  type  pattern  as  shown  for  oxygen 
stretching  region.     The  end-on  model  (be  it  C      or  C    symmetry)  is  expected 
to  have  two  closely  spaced  doublets  with  no  absorption  in  between.  The 
spacing  of  the  doublets    is  expected  to  be  small  since  only  a  small 
frequency  difference  is  to  be  expected  in  the  MO  stretching  between  the 
pairs  of  molecules  M160160  and  M*60180,  and  M180180  and  Ml80160. 

The  observed  spectrum  in  the  oxygen  stretching  region  has  been  shown 
in  Figure  1.     The  three  lines  observed  are  in  the  proper  intensity  ratios 
for  a  side-on  M02  molecule.     The  reduction  in  the  observed  frequency  from 
the  1060-1120  cm-1  region  expected  for  02  seems  to  indicate  some  charge 
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transfer  from  0=0  to  the  MO  bond.     The  narrowness  of  the  Fe160180  feature 

which  appears  to  have  a  Avi/2  of  about  1  cm-1  would  make  the  presence  of 

another  unresolved  component  to  the  doublet  (as  would  be  required  for  end- 
on  bonding)  improbable. 

The  M-0  stretching  modes  in  the  500  cm-1  region  have  already  been 
discussed.     The  weakness  of  these  modes  made  diffusion  experiments  diffi- 
cult for  all  but  the  1602  case.     It  was  also  not  possible  to  observe  all 
the  modes  of  the  three  isotopic  species  Fe1602,  Fe160180,  and  Fe1802  in 
one  experiment.     Three  separate  experiments  were  necessary  to  obtain  the 
modes  of  the  three  isotopic  species.     Despite  these  difficulties,  however 
it  should  be  noted  that  a  synthesis  of  the  three  observed  spectra  gives  a 
three  line  spectrum  with  features  at  517.1,  508.1  and  494.0  cm-1.  This 
pattern  is  expected  for  a  side-on  model  of  Fe02 •     It  should  also  be  noted 
that  the  relative  intensities  of  the  517.1,  508.1  and  494.0  cm-1  to  their 
counterparts  at  945.9,  930.8  and  911.4  cm-1  are  approximately  the  same  for 
each  pair.     This  coupled  with  the  controlled  diffusion  experiment,  possible 
with  only  the  Fe160o  sample,  showing  similar  ratios  of  disappearance  of  the 
945.9  and  517.1  cm-*  features  lends  some  confidence  to  the  assignment  of  the 
0-0  stretching  mode  and  an  Fe-0  stretching  mode  of  a  side-on  Fe02  molecule. 


A  detailed  assignment  of  the  observed  vibrations  of  Fe02  to  normal  modes 
is  not  possible  since  only  two  of  the  three  vibrations  have  been  observed. 
Prolonged  periods  of  deposition  of  Fe(g)  +  02/Ar  have  not  yielded  any  other 
features  assignable  to  Fe02 .     It  is  expected  that  the  missing  fundamental 
which  is  a  Fe-0  stretching  mode  is  at  a  lower  frequency  than  the  observed 
517.1  cm-1  for  Fe1602«     The  observation  of  only  two  vibrational  modes  is 
not  unexpected.     Of  the  previous  molecular  complexes  of  this  type  observed, 
three  modes  have  been  observed  for  only  alkali  dioxides  such  as  Li02.  Indeed 

for  most  of  the  molecules  studied  such  as  Pd09,  Ni09  etc.,   only  two  modes 
have  been  observed  for  the  M^C^  species.    Two  modes  have  also  been 
observed  for  Ba^02  and  Ba^fy^  and  only  a 

BaO  stretching  mode  for  Ba180*ko.     The  ratio  of  the  observed  frequencies  for 
the  0-0  stretching  mode  of  Fe1802  and  Fe1602  911.4/945.9  =  0.9635  is  consid- 
erably greater  than  0.9428  expected  for  a  pure  0-0  stretching  mode.  This 
indicates  that  this  vibration  is  mixing  with  an  FeO  stretching  mode.  The 
product  rule  for        v2  of  a  C^v  XY2  end-on  model  is 


i  i 

v2 
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2  my  +  m 
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The  observed  ratio  for  (v^  v2)  Fe    O2/ (v^  V2)  Fe    O2  is 


vj  v2 

~~ i  i 

vj  v2 


(911.4)  (494.0) 
(945.9)  (517.1) 


=  0.9205 
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This  is  larger  by  about  1.3%  than  the  ratio  expected  for  harmonic  vibra- 
tions.    An  alternative  assignment  of  the  494.0  and  517.1  cm-1  features 
as  due  to  V3  the  asymetric  stretching  mode  of  a  C      model  is  not  very 
satisfying.     Using  these  two  frequencies  an  O-Fe-0  angle  of  87.3°  is 
computed.     This  angle  is  larger  than  one  would  expect  for  such  a  molecule. 
A  summary  of  the  results  obtained  is  given  in  Table  1. 

There  have  been  several  solid  and  liquid  phase  studies  of  the  molecu- 
lar complexes  of  oxygen  and  hemoglobin  type  molecules.     These  studies  were 
designed  to  elucidate  the  bonding  between  oxygen  and  iron  in  hemoglobin. 
A  preliminary  report  on  an  X-ray  diffraction  study  of  a  model  porphryn 
molecule  seems  to  indicate  end-on  bonding  with  a  bent  Fe-0-0  bond  (21). 
Some  infrared  (22)  and  resonance  Raman  spectroscopic  studies  (23)  would 
seem  to  support  this  view.     Both  the  infrared  and  Raman  studies  only  inves 
tigated  1602  and  18C>2  samples.     As  was  seen  in  the  0-0  stretching  region 
previously,  both  end-on  and  side-on  bonding  WOuld  give  similar  spectra. 
The  spectra  of  the  mixed  isotope,   1&0180  bonded  to  Fe  would  clearly  dis- 
tinguish between  the  two  proposed  models  in  these  condensed  phases. 

Conclusion 

Spectra  of  the  molecular  complex  formed  by  the  reaction  Fe(g) 

with  oxygen  have  been  studied  by  the  methods  of  infrared  matrix  isolation 

spectroscopy.     Spectra  were  observed  using  both  normal  and  isotopically 

enriched  oxygen  samples.     While  a  complete  spectrum  was  not  obtainable 

due  to  the  weakness  of  the  observed  bands  the  weight  of  evidence  is 

indicative  of  a  side-on  isoceles  triangle  (C„  )  model  for  FeOo. 

zv 

Table  1 

Observed  Fundamental  Modes  for  Fe02(cm-1) 

Assignment  Fe16p?  Fe16Q?  Fe160180 

911.4  930.8 
494.0  508.1 


0-0  Stretch  945.9 
Fe-0  Stretch  517.1 
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Chapter  8 

Studies  of  the  Reactions  of  Ba(g)  with  O3 

W.  Braun,  M.  Kurylo,   S.  Abramowitz 


The  kinetics  of  the  reaction  of  Ba(g)  with  ozone  is  being  studied 
using  modulation  techniques.     By  exciting  ozone  vibrationa 1 ly  with  a 
mechanically  chopped  monochromatic  carbon  dioxide  laser,  one  can  measure 
several  parameters.     In  addition  to  being  able  to  measure  the  ratio  of 
the  rates,   k,  and  k„ 


Ba(g)  +  0^  -»    BaO  +  09  +  hv 


Ba(g)  +  03     -*    BaO  +  02  +  hv 

one  can  measure  (by  using  a  suitable  multichannel  analyzer)   the  rates 
of  growth  and  decay  of  luminescence  (produced  by  the  reaction)  in  the 
millisecond  time  regime.     In  this  way,   one  can  in  principle  measure  the. 
rates  of  production  and  loss  of  the  various  excited  species  that  may 
be  present  in  the  flame. 


The  projected  experiments  will  measure  these  variables  as  a  function 
of  argon  pressure,   in  order  to  try  to  elucidate  the  proposed  mechanisms. 
The  luminescence  will  also  be  passed  through  a  fast  monochromator  in 
order  to  investigate  its  wavelength  dependence  as  a  function  of  argon 
pressure  and  time. 

Very  preliminary  data  using  the  total  luminescence  obtainable  from 
the  reaction  indicate  a  rate  enhancement  when  vibrationally  excited  ozone 
is  used .     In  these  experiments  pressures  of  barium  in  the  ten  micron  range 
are  used.     The  barium  is  driven  by  Ar  through  a  resistively  heated 
1/4"  O.D.  stainless  steel  tube  into  a  small   reaction  chamber.     The  ozone 
which  is  generated,   using  a  commercial  ozonator,   is  vibrationally  pumped 
with  the  P.     line  of  a  CC>2   laser.     Chopping  speeds  are  used  such  that 
a  full  on-orf  wave  of  10  milliseconds  is  investigated.     Argon  pressures 
of  about  5  torr  were  utilized  in  the  initial  experiments. 
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Chapter  9 

DISSOCIATION  ENERGIES  OF  THE  SCANDIUM  GROUP  AND 
RARE-EARTH  GASEOUS  MONOXIDES 


Ralph  F.  Krause,  Jr. 


Abstract 

The  dissociation  energies  of  the  scandium  group  and 
rare-earth  gaseous  monoxides  are  evaluated  by  reviewing 
the  literature  available  on  Knudsen  effusion  rates  and 
mass-spectrometric  data.     The  rate  of  effusion  data  for 
vaporizing  the  solid  sesquioxides  to  form  gaseous  monoxides 
and  monatomic  oxygen  are  corrected  for  the  presence  of  other 
vapor  species.     These  include  metal  atoms,  dioxides,  diatomic 
oxygen,  and  oxides  of  the  container  and  except  for  the  metal 
atoms  were  not  considered  in  the  original  sources.  The 
results  from  various  sources  are  compared  by  using  the  same 
auxiliary  data  which  are  needed  to  derive  the  dissociation 
energy  of  each  monoxide.     Selected  values  of  thermodynamic 
properties  are  provided  for  all  of  the  vapor  species  and 
solid  sesquioxides  under  consideration.     Corrections  of  the 
rate  of  effusion  data  for  the  presence  of  neodymium  dioxide 
are  found  ambiguous.     Using  the  dioxide  data  from  one  source 
leads  to  a  molar  ratio  of  dioxide  to  monoxide  near  unity. 
From  another  source  the  dioxide  data,  though  inconsistent  with 
the  rate  of  effusion  data,  indicates  this  molar  ratio  to  be  negli- 
gible.    This  ambiguity  suggests  that  the  extent  of  forming  this 
dioxide  and  perhaps  others  should  be  examined  further.  When 
a  tantalum  container  is  used,  corrections  of  the  rate  of  effusion 
data  for  the  presence  of  tantalum  oxides  are  found  substantial 
in  all  the  cases  of  the  sesquioxides  effused.    When  a  tungsten 
container  is  used,  corresponding  corrections  appear  appreciable 
only  in  the  case  of  europium  sesquioxide  whose  effusion  also 
appears  to  require  an  appreciable  correction  for  diatomic  oxygen. 
Corrections  for  the  presence  of  metal  atoms  seem  negligible 
in  most  cases,  but  they  are  appreciable  in  the  effusion  of  the 
sesquioxides  of  europium,  thulium,  and  ytterbium.  Finally, 
differences  between  dissociation  energies  of  several  pairs  of 
monoxides  are  compiled  from  published  mass-spectrometric 
studies  of  isomolecular  oxygen  exchange  reactions. 
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I. 


Introduction 


The  scandium  group  and  rare-earth  metals  L  appear  as  possible 
candidates  for  chemical  lasers  in  the  visible  spectrum.  Consider 
the  bimolecular  gaseous  reaction, 


in  which  XO  may  be  N2O  or  0-,.     A  test  for  a  metal  atom  to  be  a 
candidate  is  that  the  monoxide  is  formed  in  an  excited  state  which 
has  an  excess  energy  equivalent  to  an  electronic  transition  in  the 
visible  spectrum,  namely  more  than  40  kcal  mol~^.     These  metal  monoxides 
meet  this  criterion  because  their  dissociation  energies  are  greater 
than  the  sum  of  the  above  minimum  transition  energy  and  the  dissoci- 
ation energy  of  either  XO. 

Values  of  the  dissociation  energies  D° ,  corresponding  to  the 
all  gaseous  reaction, 


are  provided  in  table  1.     The  last  column  gives  values  obtained  from 
the  present  work;  details  are  given  in  subsequent  sections.  Other 
columns  give  values  reported  earlier  in  other  compilations.  Those 
of  Suchard  (1974)  and  Rosen  (19  70),  are  both  primarily  collections 
of  spectroscopic  data.     Although  Schumm  et_  al_.   (1973)  made  selections 
from  among  published  values  which  they  adjusted  to  be  consistent 
with  their  values  of  the  heats  of  formation  of  related  compounds,  the 
documentation  of  their  efforts  has  yet  to  be  published.     Brewer  and 
Rosenblatt  (1969)  estimated  Gibbs  free  energy  functions  for  the 
gaseous  monoxides  and  recalculated  dissociation  energies  from  the 
published  values  of  selected  equilibria  measurements.     Ames  et  al. 
(1967)  performed  an  extensive  series  of  measurements.     They  used  the 
Knudsen  method  to  study  the  vaporization  of  sesquioxides  and  mass- 
spectrometry  to  study  several  isomolecular  exchange  reactions. 

Determining  accurate  dissociation  energies  from  weight-loss  and 
mass-spectrometric  data  is  fraught  by  several  difficulties.  An 
erroneous  assessment  of  the  vapor  composition  could  made  an  otherwise 
accurate  set  of  weight-loss  data  become  misleading.     It  appears  that 
earlier  investigators  failed  to  consider  adequately  the  possible 
effects  that  the  presence  of  gaseous  dioxides  might  extend.  Other 
errors  in  weight-loss  data  will  result  from  vapor  unsaturation  and 
also  from  reactions  between  the  sample  and  its  container.  Reliable 
results  from  ion  intensity  data  depend  upon  identifying  the  ion  which 
was  produced  by  electron  impact  and  properly  selecting  its  molecular 
progenitor.     While  temperature  dependent  errors  limit  the  accuracy  of 
results  obtained  by  the  2nd  Law  method,  errors  in  the  Gibbs  free  energy 
functions  will  cause  a  systematic  dispersion  of  results  obtained  by 


L  +  XO  ->■  LO    +  X 


(1) 


LO  -+  L  +  0 
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the  3rd  Law  method.     Although  It  is  beyond  the  scope  of  this  paper 
to  examine  the  extent  of  all  these  errors,  an  attempt  has  been  made 
to  become  critical  of  some. 


II.     Weight-Loss  Vaporization  of  Sesquioxides 


Several  investigators  have  studied  the  vaporization  of  the  sesqui- 
oxides by  weight-loss  methods,  including  Knudsen-ef fusion  in  tungsten, 
tantalum,  or  rhenium  cells,  Langrauir-tors ion  in  molybdenum  cages,  and 
entrainment  in  a  container  of  the  sample  itself.     In  many  cases  inves- 
tigators had  assumed  that  the  vaporization  process  occurred  practically 
as 


Even  though  reaction  3  may  not  be  a  practical  process  for  evaluating 
the  weight-loss  data  of  all  the  sesquioxides,  it  will  be  taken  in  this 
paper  as  a  convenient  reaction  whose  equilibrium  will  be  evaluated. 

Vapor  pressure  measurements  as  a  function  of  temperature  are 
meaningful  if  the  chemical  system  Is  univariant.     When  the  sample  container 
is  excluded,  one  may  realize  from  the  Gibbs  Phase  Rule  that  a  two  com- 
ponent system  will  be  univariant  if  three  phases,  for  example  two  solids 
and  the  gas  phase,  are  present.    Although  White  et  al.  (1962)  perceived 
by  x-ray  powder  analysis  only  a  single  solid  phase  during  the  complete 
process  of  vaporizing  each  of  several  sesquioxides,  the  composition  of 
this  phase  remained  constant  for  a  particular  sesquioxide    and  varied 
among  them  from  2.96  to  2.99  atomic  oxygen  per  2  L.     White  et  al. 
also  observed  constant  effusion  rates  between  the  limits  of  20  and 
80%  of  sample  vaporized.     Consequently,  it  appears  that  these  sesqui- 
oxides vaporize  as  azeotropes,  making  a  two  component  system  univariant. 
If  the  solid-phase  were  a  mixture  of  L2O3  and  LO,  the  mol  fraction  of 
the  latter  would  appear  as  high  as  0.077;  and  if  the  solid  phase  were  an 
ideal  solution,  this  mol  fraction  would  be  the  error  in  the  pressure  of 
L0(g).     It  can  be  shown  that  for  measurements  at  2500  K  a  correction 
for  this  error  would  lower  the  calculated  dissociation  energy  by  0.6 
kcal  mol~l. 

In  accord  with  reaction  3,  weight-loss  data  alone  would  uniquely 
define  the  vapor  composition  if  the  vapor  contained  only  two  gas 
species.     In  general, "the  vapor  may  be  considered  to  contain  many 
species -for  which  an  investigator  should  account.     When  a  sesquioxide 
is,  for  example,  in  a  tungsten  container,  vaporization  may  be  expressed 
as 


L203(s)  ->  2  L0(g)  +  0(g). 


(3) 


L203(s)  +  n  WOi)  ">  a  0(g)  +  b  02(g)  +  c  L(g)  +  d  L0(g)  + 
e  L02(g)  +  T.  k±  Wx0  (g)  +  ... 


(4) 


If  reaction  A  is  assumed  to  vaporize  congruently,  a  mass  balance 
provides 


c  +  d  +  e  + 


=  2 


(5) 
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a  +  2b-fd  +  2e  +  £  y  k    +  . . .  -  3 


(6) 


and 

£  xjLk1  =  n  (7) 

The  vapor  composition  of  more  than  two  gas  species  is  not  defined  by 
weight-loss  data  unless  other  data  are  available  to  relate  all  the 
variables.     Table  2  gives  thermodynamic  values  for  several  equilibria 
which  involve  gas  species  implicated  by  reaction  4. 

One  of  the  gas  species  which  remains  undefined  is  the  metal  atom. 
Table  3  shows  that  the  mass-spectrometric  ion-intensities  from  sesqui- 
oxide  vaporization  indicate  that  the  L  /L0+  ratio  has  two  cycles  of 
increasing  values  for  the  series  of  rare  earths  including  lanthanum. 
These  values  might  be  on  the  high  side  because  the  high  ionization 
potentials  which  were  used  might  have  fragmented  the  L0+.  Ackerraann 
et  al.   (1964)  point  out  that  some  caution  must  be  exercised  in  evalu- 
ating a  meaningful  L+/LO+  ratio,  which  they  believe  cannot  be  obtained 
at  a  single  arbitrarily  selected  electron  energy.     Using  ionization 
efficiency  curves,  they  observed  a  saturation  plateau  for  Y+  at  approxi 
mately  10-12  eV  and  for  Y0+  at  approximately  20  eV;  from  these  they 
evaluated  the  Y+/Y0+  saturation  ratio    to  be  0.015,  which  is  about  1/4 
that  observed  by  Ames  e_t  al.     Extrapolating  these  curves  to  zero 
intensity  they  found  appearance  potentials  of  5.5  and  6.0  eV  for  Y0+ 
and  Y+,  respectively.     If  they  had  observed  an  appearance  potential  for 
Y      near  13  eV,  one  would  have  suspected  its  presence  due  solely  to 
fragmentation  of  Y0+,  which  their  ionization  efficiency  curve  showed 
to  be  appreciable  above  20  eV.     Conversion  of  a  ratio  of  ion-intensitie 
to  a  ratio  of  pressures  in  the  last  column  of  table  3  depends  on  the 
proportionality  of  the  following  relation 

PjL  «  I1+T/ai  (8) 

in  which        is  a  cross-section  for  detection  of  the  species  i  in  a 
particular  instrument.  +Using  rate  of  effusion  data,  Ames  et  al. 
(1967)  deduced        for  Y  ,  YO  ,  La+.  and  LaO+  with  respect  to  that 
of  Ag+.     The  results  for  the  L+/L0    ratio  of  cross-sections  averaged 
1.6  with  a  difference  of  0.2.     Although  the  values  in  the  last  column 
of  table  3  may  be  useful  for  comparison,  a  preferrable  method,  which 
is  similar  to  one  used  by  White  et  al .   (1962),  provides  a  relationship 
between  the  pressure  of  a  metal  atom  and  the  published  data  for  the 
atomization  of  its  sesquioxide. 

L203(s)  +  2  L(g)  +  3  0(g)  (9) 

and  this  method  will  be  described  below. 

The  weight-loss  method  which  was  predominantly  used  for  the  sesqui 
oxide  vaporization  measurements  in  the  literature  was  the  Knudsen 
effusion  method.    At  temperature  _T  the  vapor  pressure        for  a  vapor 
species  of  molecular  weight        is  obtained  from 
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Table  2*     Some  Equilibria  Involved  with  the  Vaporization  of  Sesquioxides. 


Reaction6  AHq   A[-(G°-Hn)/Tl/cal  moj~V^  Ref« 


Vpal     TTinl  ^" 
i\  La  1  UlUl 

2000  K 

2200  K 

2400  K 

2600  K 

2  0  ■+  02 

-117.98 

-29.91 

-30.10 

-30.26 

OA  /.I 

-30. 41 

a 

Nd  +  Nd02  ■*  2  NdO 

-  10.5 

4.01 

3.96 

3.65 

3.07 

b 

NdC>2  ■*■  NdO  +  0 

c 

101 

26.93 

27.01 

27 .06 

Z/.ll 

Ta(s)  +  0  +  TaO 

-  6.88 

7.24 

6.95 

6.69 

6.44 

d 

Ta(s)  +20+  Ta02 

-163.93 

-27.10 

-27.41 

-27.67 

-27.90 

d 

W(s)  +  0  -y  WO 

42.74 

9.41 

9.11 

8.84 

8.57 

a 

W(s)  +  2  0  ->  WO 

L 

-99.10 

-20.02 

-20.41 

-20.76 

-21.08 

a 

W(s)  +  3  0  -*  W03 

-245.92 

-57.65 

-58.04 

-58.38 

-58.69 

a 

2  W(s)  +  60  +  (WO.). 

-630.70 

-146.71 

-147.03 

-147.31 

-147.56 

a 

3  W(s)  +80+  W  Og 

-878.30 

-212.72 

-213.15 

-213.51 

-213.84 

a 

3  W(s)  +90+  (W03)3 

-1011.56 

-247.87 

-248.30 

-248.65 

-248.97 

a 

4  W(s)  +  12  0  +  (WOa). 

-1374.59 

-344.85 

-345.27 

-345.61 

-345.92 

a 

aStull  and  Prophet  (1971). 
Pupp  and  Gingerich  (1971). 

Calculated  by  the  3rd  Law  from  the  Kp  data  of  Staley  and  Norman  (1969) . 
dSchick  (1966). 
Gaseous  species  except  where  noted. 


Table  3.     Ratio  of  Mass-Spectrometric  Ion  Intensities  from  Vaporization  of  the  Sesquioxides. 

T  


Solid 


Panish  (1961),       Shchukarev  &  Semenov  (19617 
20  eV  10  eV  45  eV 


Ames  et.  al.   (196  Tj~ 
25  eV 


Average 


PT/P  a 
L  LO 


SC2°3 

l/8b 

0.2 

0.16 

0.1 

Y2°3 

1/16 

0.06 

0.04 

La2°3 

0.16 

<0.01 

0.08 

0.05 

Ce2°3 

0.20 

0.2 

0.12 

Pr203 

<0.1 

0.36 

0.2 

0.1 

Nd203 

<0.1 

0.25 

<0.01 

0.1 

0.06 

Sm203 

0.5  to  1.0 

1.3 

1.0 

0.5 

0.9 

0.6 

EU2°3 

8±2,  9±1 

7 

5 

10 

9 

6 

Gd203 

<0.1 

1.2 

0.3 

0.5 

0.3 

Tb2°3 

<0.1 

0.9 

0.2 

0.4 

0.3 

Dy2o3 

0.4 

1.7 

0.5 

0.5 

0.8 

0.5 

H°2°3 

1.0 

2.2 

1.3 

0.5 

1.2 

0.8 

Er203 

1.5  ±  0.5 

1.4 

1.0 

0.5 

1.1 

0.7 

Tm203 

15  ±  5 

10 

10 

12 

8 

Yb203 

(>30) 

>50 

50 

30 

Lu203  . 

<0.1 

1.7 

1.1 

0.5 

0.8 

0.5 

The  ratio  of  cross-sections 

(L+/L0+) 

was  assumed  as  1.6; 

see  text. 

^Semenov 

(1965). 
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P    =  wiV2  RT/M  /(tBWbf)  (10) 

in  which        is  the  mass  of  species  i  that  effuses  in  time  _t  through 
an  orifice  with  a  Clausing  factor        and  an  area  IS.     The  value  of 
the  constant  R  is  (8.  314  34)107  erg  K-lmol-1,  and  1  atra  is  (1.01325)106 
dyn  cm~2.     Paule  and  Margrave  (1967)  show  that  the  fraction  of 
saturation  _f  of  the  vapor  in  a  Knudsen  cell  with  a  Clausing  factor 
W_  and  an  area  A  may  be  estimated  by 

f  =  [1  +  (BW,/A)(a_1  +  W"1  -2)1_1  (11) 
D  a 

in  which  _a  is  the  vaporization  coefficient  of  the  condensed  sample. 

Values  of  Clausing  factors  are  tabulated  in  an  appendix  of  reference 

35  as  functions  of  the  length  and  radius  of  a  cylinder.     The  mass  w^ 

of  species  i^  is  related  to  the  total  mass  w  of  sesquioxide  of  molecular 
weight  M  by 

w±  -  w  k  M1/M  (12) 

in  which  k^  is  the  coefficient  of  a  vapor  species  in  reaction  4. 
When  eq.  10  is  used  to  calculate  a  hypothetical  Knudsen  corres- 
ponding to  the  w  and  M  of  a  sesquioxide, 


Pi  =  ki  PkVVM  (13) 

Using  eq.  13,  the  k^  of  species  i  may  be  expressed  as  a  function  of 
a  known  Kp  for  a  reaction  in  which  the  species  i  is  also  involved;  for 
example  in  the  case  of  a  metal  atom  L  in  reaction  9 


c    =  Kp  ^    /(a  Pk  ^      )  (14) 


The  K    of  a  reaction  at  a  given  temperature  is  calculated  from  the 
corresponding  AHg  and  A [-(G°-Hq) /T]  as 

K    =  exp{A[-(G°-H°)/T]/R  -  AH«/(RT)}  (15) 
p  u  u 

In  like  manner  the  procedure  of  the  preceding  example  may  be  applied 
to  other  equilibria,  such  as  those  in  table  2,  which  contribute  to 
the  definition  of  vapor  species  in  reaction  4.    After  values  of  a  and 
<i  are  evaluated  by  simultaneous  solution  of  eqs.  5  and  6,  the  Kp  of 
reaction  3  may  be  calculated  simply  as 

Kp(3)  '  PL02  P0  (16) 
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in  which  the  _Pjl  are  determined  by  eq.  13.     When  the  Langtnuir  method 
is  employed,  eq.  13  could  also  be  used  since  it  involves  using  the 
Y^i  term  of  eq.  10  similar  to  the  Knudsen  method.     If  reliable  Ffc 
data  as  a  function  of  temperature  are  available,  the  d(£n  P]c)/d(l/T) 
could  be  applied  to  define  an  unknown  variable  in  eqs.  5  and  6. 

In  principle  the  entrainment  method  offers  yet  another  means  to 
evaluate  an  unknown  variable  in  eqs.  5  and  6  when  it  is  combined  with 
the  Knudsen  method.     In  practice,  however,  experimental  imprecision  of 
the  two  methods  often  outweighs  any  meaningful  contribution.  The 
entrainment  method,  which  is  truly  a  vapor  density  measurement,  would 
result  in  defining  the  pressure  of  vapor  species  i_  as 

Pt  =  (w  /V)RT/M  (17) 

When  eq.   17  is  used  to  calculate  a  hypothetical  entrainment  Pg 
corresponding  to  the  w  and  M  of  a  sesquioxide, 

Pt  =  kt  Pe  (18) 

which  is  analogous  to  eq.  13. 

III.    Auxiliary  Data 


The  published  data  necessary  to  evaluate  the  Kp  of  reaction  9  are 
assembled  in  tables  4  and  5.     Holley  et  al.   (1968)  nad  reviewed 
extensively  the  published  thermodynamic  data  for  the  sesquioxides. 
Table  4  includes  some  revised  values  which  have  been  measured  recently. 
The  values  of  (H°-Hq)  and  -(G0-Hq)/T  in  table  4  are  extrapolations 
which  were  calculated  using  expressions  of  the  form 

(H°-H°     )  =  A  +  BT  +  CT2  +  D/T  (19) 

which  Holley  et  a_l.  had  determined  from  relative  enthalpy  data. 
Holley  e_t  al^.  had  selected  the  coefficients  J)  so  that  the  high  temp- 
erature data  joined  smoothly  with  the  low  temperature  heat  capacity 
data.     Since  the  sources  in  table  5  gave  values  of  the  Gibbs  free 
energy  functions  with  respect  to  298.15  K,   they  were  appropriately 
adjusted  to  zero  Kelvin. 

Since  reaction  9  is  the  sum  of  reaction  3  and  twice  reaction  2, 
the  dissociation  energy  of  a  monoxide  is  related  as 

D°(L0)  =  [AH°(9)  -  AH°(3)]/2  (20) 
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The  2nd  or  3rd  Law  Method  may  be  used  to  obtain  the  of  reaction  3 

from  measured  values  of  the  corresponding^.     I  f  _Kp  data  are  available 
as  a  function  of  temperature,   the  2nd  Law  would  give 

AH°  -  -ft[d(in  Kp)/d(l/m  -  A.(Hj-H°)  (21) 

in  which  values  of  (H™-H°)   for  atomic  oxygen  were  obtained  from  Stull 
and  Prophet  (1970)  ana  values  of  (H£-H°)   for  all  these  monoxides  were 
obtained  from  Ames  et.  a_l.  as  16.82,    18.59,   20.37,   and  22.15  cal  mol"LK 
at  2000,  2200,  2400,  and  2600  K,   respectively.     At  a  given  temperature 
the  3rd  Law  gives 

AHj  =  T  |a[-(G°-H°)/T]-  R  in  K  }  (22) 

in  which  values  of  -(G°-Hq)/T  for  all  these  monoxides  were  also 
obtained  from  Ames  e_t.  a_l.  and  are  given  in  table  6.     Differences  in 
the  results  of  the  two  methods  indicate  either  errors  in  K  ,   in  its 
derivative  with  temperature,  or  in  the  auxiliary  thermodynamic  functions. 

The  -(G°-Hq)/T  values  in  table  6  are  those  used  in  the  present  work, 
but  they  appear    somewhat  uncertain.     Table  7  shows  -(G°-H°)/T  values 
which  were  estimated  either  by  Ames  e_t.  al .   (1967)  or  by  Brewer  and 
Rosenblatt  (1969).     Using  the  values  of  the  latter  would  lead  to  values 
of  dissociation  energies  which  are  approximately  6  kcal  mol"  smaller 
than  those  derived  from  table  6.     Both  sources  appear  to  have  selected, 
respectively,   similar  ground  state  vibrational  frequencies  for  all  the 
rare-earth  monoxides,  and  these  estimates  are  approximately  the  same  as 
those  which  DeKock  and  Weltner  (1971)  observed  recently  in  the  infrared 
between  808  and  832  cm"    by  matrix  isolation. 

The  major  cause  of  the  discrepancy  is  that  the  above  two  sources 
chose  different  electronic  multiplicities  for  the  ground  and  low- 
lying  states.     These  multiplicities  affect  mainly  the  entropy  of  a 
monoxide.    While  both  sources  reported  using  a  divalent-ion  model  for 
the  electronic  partition  functions,  Brewer  and  Rosenblatt  used  higher 
multiplicity  values;   for  example,  with  ScO  Ames  e_t.a_l.  considered  only 
the  ground  state  and  assumed  it  to  have  a  multiplicity  of  2,  whereas 
Brewer  et .  al .  used  4  along  with  that  of  6  for  an  assumed  low-lying 
energy  level  which  they  took  as  197  cm"*,   the  value  for  the  divalent 
ion.     The  values  of  Ames  et .  a  1 .  are  nevertheless  preferred.  Approximate 
wavefunctions  for  ScO  were  calculated  by  Carlson  e_t.  al .   (1965)  who 
concluded  that  2^  is  the  ground  state  of  the  molecule.     Kasai  and 
Weltner  (1965)  presented  experimental  evidence  which  they  believed  was 
conclusive  that  ScO,  YO,   and  LaO  all  have        ground  states.  Finally 
Green  (1971)  presented  evidence  that  the  unobserved  low-lying  electronic 
states  of  ScO,  YO,   and  LaO  can  be  interpreted  using  a  model,   but  he 
concluded  that  these  states  make  no  appreciable  contribution  to  the 
electronic  partition  functions. 
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Table  6. 

Values  of  the 
Monoxides . 

Gibbs  Free  Energy  Functions 

for  the 

Gaseous 

Oxide 

Electronic 

-(G*-H°)/t/  cal  Tnol_1K 

-1 

Multiplicity 

2000  K 

2200  K 

2400  K 

2600  K 

ScO 

2 

61.19 

61.99 

.  62. 72 

A 1  10 

YO 

2 

63.55 

64.  35 

65 .09 

A  ^  77 

LaO 

2 

64.97 

65 . 78 

66. 51 

£7  on 
o  /  •  zu 

CeO 

3 

65.71 

66.52 

67.25 

0/i7J 

PrO 

4 

66.37 

67. 18 

67.92 

Aft  AH 

NdO 

5 

66.86 

67.67 

68.41 

SmO 

7 

67.66 

68.46 

69.20 

AQ    ft  A 

EuO 

8 

67.96 

68.  76 

69.50 

/ U  •  lu 

GdO 

9 

68.29 

69.09 

69.83 

70  M 

TbO 

8 

68.09 

68.90 

69.63 

7fi  11 

DyO 

7 

67.89 

68.  70 

69 .4  3 

HoO 

6 

67.63 

68.43 

69.17 

AQ  ft^ 

ErO 

5 

67.31 

68.11 

68.85 

AQ    ^  1 

TmO 

4 

66.90 

67.70 

68.44 

AQ    1  9 

YbO 

3 

66.40 

A  7  on 

0  /  .  74 

Aft    A  9 

LuO 

2 

65.61 

66.42 

67.15 

67.83 

Reference 

.    Ame  s  e  t .  a  1 

.  (1967). 

Table  7.     Differences  in  Dissociation  Energies  from  Different  Values  of  the 
Gibbs  Free  Energy  Functions  for  the  Gaseous  Monoxides. 


Oxide 

T/K 

-<G°-H°)/y 

cal  mol  K 

AD°(L0)a 

0 

Ames  et.  al.  (1967) 

Brewer  et.  al.  (1969) 

kcal  mol  1 

ScO 

2600 

63.39 

66.46 

8.0 

YO 

2600 

65.77 

68.53 

7.2 

LaO 

2400 

66.51 

68.79 

..  5.5 

CeO 

2000 

65.71 

68.90 

6.4 

PrO 

2000 

66.  37 

53. 0b 

NdO 

2400 

68.41 

54. 4b 

SmO 

2400 

69.20 

71.5 

5.5 

EuO 

2200 

68.76 

70.7 

4.3 

GdO 

2600 

70.51 

72.55 

5.3 

TbO 

2000 

68.09 

69.2 

2.2 

DyO 

2600 

70.11 

72.4 

6.0 

HoO 

2600 

69.85 

72.2 

6.1 

ErO 

2600 

69.53 

72.0 

6.4 

TmO 

2400 

68.44 

71.0 

6.1 

YbO 

2400 

67.94 

70.6 

6.4 

LuO 

2600 

67.83 

53. 6b 

'Calculated  as  Ames'  value  of  [-(G°-H°)/T]  less  Brewer's  value. 

o 

Apparently  an  error  in  the  source's  text. 

134 


IV.     Errors  in  Vapor  Composition 


A  possibly  serious  systematic  error  in  the  calculations  using 
weight-loss  data  might  exist  if  the  formation  of  metal  dioxides  and 
more  complex  oxides  are  neglected.     As  shown  above  in  table  2,  Staley 
and  Norman  (1969)  and  Pupp  and  Gingerich  (1971)  have  studied  equilibria 
involving  NdC^.     Since  their  data  will  be  shown  to  be  disparate,  the 
extent  of  Nd02  formation  is  difficult  to  evaluate.     The  last  column  of 
table  8  shows  values  of  Dq  for  NdO  which  were  calculated  by  varying 
the  AHq  as  shown  in  column  1.     Corrections  were  applied  for  the  formation 
of  diatomic  oxygen,  tungsten  oxides,  and  the  Nd  atom,  using  the 
sesquioxide  atomization  data.     A  value  of  AHq  for  reaction  c_  that  is 
as  low  as  that  shown  in  the  last  line  of  column  2  is  inconsistent 
with  the  rate  of  effusion  data  because  there  is  produced  P^d/^NdO  < 
0.     In  the  next  to  the  last  line  the  mo 1  fraction  of  Nd02  is  assumed  zero. 
If  the  data  of  Pupp  and  Gingerich  were  believed,  one  would  judge  from 
the  results  in  line  1  that  Nd203  vaporizes  practically  as 

Nd203(s)  ■>  NdO(g)  +  Nd02(g)  (23) 

Also,  the  negligible  amount  of  atomic  oxygen  shown  in  line  1  would 
correspond  to  the  fact  that  no  tungsten  oxides  were  detected  mass- 
spectrometrically  by  Goldstein  e^  a_l.   (1961)  .     The  presence  of  other 
metal  dioxides  are  apparent.     De  Kock  and  Weltner  have  observed 
infrared  spectra  which  they  have  assigned  to  CeCb,  Pr(^2»  an<^  TbC^. 
The  extent  of  formation  of  these  and  other  dioxides  would  help  define 
the  vapor  composition  of  weight-loss  studies  and  provide  more 
accurately  calculated  dissociation  energies. 

The  use  of  mass-spectrometry  to  study  equilibria  involving 
LO2  might  entail  some  difficulties  especially  when  too  high  an 
ionization  potential  is  used.     Using  ionization  efficiency  curves 
for  CeO+  and  Ce02+,  Ackermann  and  Rauh  (1*571),  ref.   3,  found  appearance 
potentials  of  5.2  and  10.3  eV,  respectively,  based  on  the  known  5.5  eV 
ionization  potential  for  Ce+.     They  found  the  onset  of  fragmentation 
of  Ce02  to  CeO+  at  about  11  eV  and  of  CeO  to  Ce+  at  about+13.5  eV. 
After  an  interval  of  time  they  observed  the  amoun£  of  CeO  had 
increased  by  a  factor  of  7  and  the  amount  of  Ce02    had  decreased  by 
a  factor  of  2.     A  decrease  in  the  pressure  of  atomic  oxygen  would 
account  for  the  shift  of  Ce02  to  CeO.     It  seems  entirely  possible 
that  the  use  of  mass-spectrometry  could  miss  LO2  entirely  unless 
a  sufficient  pressure  of  oxygen  were  present. 

Errors  which  result  from  sample  reaction  with  its  container  appear 
most  pronounced  when  sesquioxides  are  vaporized  in  tantalum  Knudsen 
cells.     Using  the  rate  of  effusion  data  of  Kulvarskaia  and  Maslovskaia 
(1960),  eqs.  5  and  6  were  solved  simultaneously  with  and  without 
corrections  for  the  tantalum  oxides.     The  results  are  shown  in  table  9. 
Corrections  for  diatomic  oxygen  and  the  sample-metal  atoms,  using  the 
atomization  data  for  the  sesquioxides,  were  applied;  but  the  sample- 
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Table  8.     Effect    of  the  Uncertainty  of  Forming  Neodymium 
Dioxide  Gas  on  Calculations  Using  Sesquioxide  Rate  of  Effusion  Datac 


Assumed 

Calculated 

Averaged  Results 

( 30  runs , 

2255-2434  K) 

AH°  (b) 

AH°(c) 

P  /P 

Mole  Ratio  per  Nd,, 

0^  (vaporized) 

D°(Nd0) f 

kcal  niol  ^ 

kcal  mo  1  ^" 

Nd  NdO 

W  (vaporized) 

0(g) 

NdO(g) 

NdO  (e) 

kcal  mol  ^ 

-iod 

148.  3 

0.016 

/„* ,„-5 
(2)10 

0.007 

0.9  7 

1.01 

158.2 

-20 

141.7 

0.002 

(9)10-5 

0.026 

1.02 

0.98 

16l!7 

-30 

135.0 

(3)10-^ 

c\  no/. 

1.09 

0.91 

165.0 

-40 

128.2 

(5)10"° 

0.003 

0.27 

1.28 

0.72 

168.2 

-50 

120.9 

(1.4)10" 

5  0.012 

0.56 

1.60 

0.40 

170.9 

(6)10-5 

0.03 

0.87 

2.00 

0 

173.0 

101e 

<0 

aGoldstein  et 

al.  (1961). 

C  Nd02(g)  ■>  NdO(g)  + 

0(g) 

SStaley  and 

Norman  (1969) 

bNd(g)  +  Nd02 

(g)  -  2  NdO(g) 

d 

Pupp  and  Gingerich 

(1971) 

fThird  Law 

values 

Table  9.     Effects  of  Correcting  for  Gaseous  Tantalum  Oxides  on  Calculations  Using  Sesquioxide  Rate  of  Effusion  Datac 


L203(s) 

Without  Correction 

With  Correction 

PL/PL0 

Mole  Ratio 
0(g)/L203(vaporized) 

D°(LO)b 

VPL0 

Mole  Ratio  per  L203 

(vaporized) 

D°(L0)b 

kcal  mol  ^ 

Ta(vaporized) 

0(g) 

kcal  mol  ^ 

La2°3 

io-9 

0.61 

208 

io"6 

0.58 

0.0033 

197 

Pr2°3 

io-6 

0.85 

190 

io-A 

0.68 

0.009 

180 

Nd2°3 

IO"5 

0.87 

179 

io"2 

0.72 

0.014 

169 

Sm^ 

io"4 

0.8'0 

154 

0.28 

0.90 

0.009 

143 

Eu2°3 

io'4 

0.41 

140 

2.0 

1.23 

0.0032 

125 

177 

Gd2°3 

IO"6 

0.91 

187 

io"3 

0.72 

0.011 

Dy203 

io"4 

0.94 

164 

0.11 

0.94 

0.021 

155 

Ho203 

io"A 

0.97 

166 

0.13 

1.03 

0.023 

156 

Er2°3 

io"3 

0.99 

166 

0.16 

1.12 

0.027 

157 

Yb203 

io"3 

0.76 

138 

3.1 

1.44 

0.010 

123 

Lu203 

io"6 

0.90 

19  3 

io-A 

0.76 

0.012 

183 

^ulvarskaia  and  Maslovskaia  (I960) . 


Third  Law  values. 
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metal  dioxides  were  assumed  absent.     The  general  effect  of  correcting 
for  the  tantalum  oxides  is  to  reduce  the  value  of  the  calculated 
dissociation  energy  by  approximately  10  kcal  mol~l.     The  correction 
of  course  depends  upon  the  pressure  of  atomic  oxygen;  so  if  the 
formation  of  a  sample-metal  dioxide  were  appreciable,  alternative 
results  would  be  obtained.     Goldstein  e_t  a_l.   (I960)  detected  no 
tantalum  oxides  mass-spectrometrically  when  they  s tud led  sesquioxide 
effusion  in  a  tantalum  cell,  but  they  observed  that  the  total  weight- 
loss  of  the  cell  and  sample  exceeded  the  sample  mass. 

The  effects  that  other  assumptions  have  upon  the  calculations  of 
dissociation  energies  are  shown  in  table  10.     Here  are  analysed  the 
sesquioxide  rate  of  the  effusion  data  of  Goldstein  e_t  al_.   (1961)  and 
of  Ames  et  al,   (1967),  both  of  whom  had  used  a  tungsten  cell.  The 
variation  in  the  calculated  dissociation  energy  for  LaO  is 
representative  of  those  for  the  monoxides  of  Sc,  Y,  Nd,  Sm,  Gd,  Dy, 
Ho,  Er,  and  Lu.     Calculations  using  the  effusion  data  for  YD2O3 
yielded  negative  P(L)/P(L0).     In  assumptions  1  through  4,   the  sample- 
metal  atoms  were  corrected  by  the  atomization  data  for  the  sesquioxides. 
The  first  assumption,  which  is  the  one  employed  by  Ames  et  al., 
considers  that  only  the  a,  £,  and  d  of  reaction  4  are  important. 
The  second  assumption  includes  _b,  and  the  third  assumption  includes 
the        for  the  various  tungsten  oxides.     Assumptions  3  and  5  differ 
in  that  the  latter  uses  values  in  the  last  column  of  table  3  instead 
of  the  atomization  data.     The  P(L)/P(L0)   calculated  by  assumption  3 
for  all  the  sesquioxides  follows  the  same  two  cvclic  pattern  as 
table  3  except  that  assumption  3  gives  them  one  or  more  orders  of 
magnitude  smaller.     Even  though  Ames  ejt  al_.  detected  no  tungsten  oxides 
mass-spectrometrically,  they  did  observe  nearby  0.1  mole  of 
tungsten-loss  per  mole  of  L2O3.     The  results  in  table  10  indicate 
that  the  corrections  for  the  tungsten  oxides  have  appreciable  effects 
mainly  for  the  case  of  E^O^.     This  correction  is  confirmed  to  some 
extent  by  the  results  of  Petzel  and  Greis  (1971)  who  used  iriduim- 
lined  Knudsen  cells  to  study  the  vaporixation  of  EU2O3  in  the  presence 
of  W. 

An  approximation  of  the  effect  of  vapor  unsaturation  in  the  Knudsen 
cell  is  shown  by  assumption  4  in  table  10.     Although  vapor  unsaturation 
produces  a  shift  in  the  vapor  equilibria  of  complex  vaporization,  the 
approximation  is  based  upon  an  estimate  of      by  eq.  11.     The  cell 
dimensions  reported  by  Goldstein  e_t  a_l.   (1961)  and  used  by  Ames  e_t  al. 
were  such  that  A/B ,  V^,  and       were  evaluated  as  16,  0.3,  and  0.5, 
respectively.     Goldstein  et  al.  (1961)  believed  that  a  >  0.5;  so  it 
appears  that  they  and  Ames  e_t  al.  reported  mass-effusion  rates  which 
included  a  value  of  _f  as  0.91.     More  recently  Messier  (1967)  observed 
Knudsen  effusion  rates  for  Gd203  through  each  of  two  orifices  of 
different  size.    As  will  be  described  later  in  Section  VI  under  Messier, 
the  vaporization  coefficient  0^  for  r^O-^  was  deduced  to  be  0.074. 
When  this  value  of  _a  is  used  to  revise  the  work  of  Goldstein  et  al., 
_f  becomes  0.68  instead  of  0.91  to  produce  the  results  of  assumption  4. 
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V.     Equilibria  Data  by  Mass  Spectrometry 


Differences  between  the  dissociation  energies  of  several  pairs 
of  monoxides,  LO  and  MO,  have  been  determined  from  equilibria  studies 
of  isomolecular  oxygen  exchange  reactions  in  the  gas  phase 

L  +  MO  ■+  LO  +  M  (24) 

Several  investigators  have  used  mass-spectrometry  to  measure  ion- 
intensities  as  a  function  of  temperature  for  different  condensed 
mixtures  of  two  or  more  metals  and  oxides,  or  of  oxides  only.  They 
calculated  the        of  reaction  24  by  assuming  that  the  corresponding 
ratio  of  cross-sections  were  unity,  which  seems  satisfactory  according 
to  the  measurements  of  the  cross-sections  for  Y+,  Y0+,  La+,  and  LaO+ 
discussed  in  Section  II  above.     Results  by  the  2nd  and  3rd  Law  methods 
are  assembled  in  table  11.     The  extent  of  fragmentation  of  LO  to  L 
was  neglected;  however,  Smoes  e_t  al_.   (1965)  and  Ackermann  and  Rauh  (ref.  2) 
were  so  concerned  about  the  ill-effects  of  fragmentation  that  they 
used  lower  ionization  potentials.     Values  of  the  dissociation  energies 
in  the  last  column  of  table  11  were  calculated  from  the  selected 
values  in  the  next  to  the  last  column.     The  results  for  LaO  has 
particular  significance  because  it  depends  upon  that  for  SiO  which 
was  derived  from  an  independent  review  by  Stull  and  Prophet  (1971). 

Another  method  which  has  been  used  to  evaluate  dissociaiton 
energies  by  mass-spectrometry  involves  measuring  the  LO VL  ratio 
of  ion- in  tens ities  for  the  vapor  species  which  are  in  thermodynamic 
equilibrium  with  both  L(£)  and  L203(s) .    A  difficulty  with  this 
method  is  the  uncertain  activities  of  the  condensed  phases,  i.e. 
their  mutual  solubilities.     Using  eq.  8  and  independent  data  for 
the  pressure  of  L(g) ,  the  Kp  of  reaction  25 

1/3  L(g)  +  1/3  L203(s)  ->  LO(g)  (25) 

may  be  calculated  from 

%  =PL2/3  (ILO+/IL+)(oi>,.0+>  <26) 

in  which  the  L+/LO+  ratio  of  cross-sections  was  shown  in  Section  II 
above  to  be  approximately  1.6. 

A  third  application  of  mass-spectrometry  involves  measuring  as 
a  function  of  temperature  the  L0+  ion-intensity  which  is  proportional 
to  the  pressure  of  LO  vapor  in  equilibrium  with  L20^(s).    When  the 
pressure  of  atomic  oxygen  is  assumed  proportional  to  the  pressure  of 
the  monoxide  over  the  temperature  interval  of  measurement,  the 
enthalpy  change  of  reaction  3  may  be  approximated  from 

AH°  =  -3R[d(ln  PL0)/d(l/T)]  (27) 
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VI.     Recalculated  Results 


Using  the  weight-loss  and  ion-intensity  data  which  are  available 
in  the  literature,  the  K    of  reactions  3  and  25  are  evaluated;  and  the 
corresponding  dissociatTEn  energies  are  calculated  from  the  auxiliary 
data  which  are  given  in  Section  III.     Second  Law  and  average  Third 
Law  values  of  the  dissociation  energies  for  the  monoxides  are  reported 
in  table  12.     Selections  of  Dq  values  listed  in  tables  11  and  12,  are 
compiled  in  the  last  column  of  table  1.     Comments  on  the  results 
from  each  source  follow  below  in  the  order  that  they  appeared  in  the 
literature. 

Chupka  e_t  aJL.   (1956)  had  used  mass-spectrometry  in  4  experiments 
to  measure  the  LaO+/La+  ion-intensity  ratio,  which  is    proportional  to 
the  LaO/La  pressure  ratio  in  equilibrium  with  La(&)  and  La20^(s) . 
Instead  of  using  their  estimate  of  0.5,  the  value  of  1.6  deduced 
by  Ames  e_t  a_l.  was  used  in  eq.  26  for  the  La+/LaO+  ratio  of  cross- 
sections.     The  resulting  Dq  value  is  exceptionally  close  to  that  in 
table  11  for  the  isomolecular  reaction  involving  Si. 

Kulvarskaia  and  Maslovskaia  (I960)  had  studied  the  vaporization 
of  several  sesquioxides  by  Knudsen  effusion.     The  3rd  Law  Dq  in  table  12 
are  those  from  table  9  which  were  corrected  for  tantalum  oxides. 
Second    Law  values  of  Dq  were  obtained  by  using  eqs.  20  and  21. 

Goldstein  e_t  al .  (1961)  had  used  both  Knudsen-ef fusion  and  mass- 
spectrometry  to  study  the  vaporization  of  La203  and  Nd203.     The  Knudsen 
value  of  the  dissociation  energy  for  NdO  was  selected  as  the  average 
of  the  extremes  which  are  reported  in  table  8  and  is  strikingly  close 
to  the  mass-spectrometric  value  which  resulted,  using  eq.  27.  The 
Knudsen  value  of  Dq  for  LaO  is  the  result  of  assumption  3  described 
above  for  table  10.     Similar  to  the  NdO  results  the  higher  value  of 
Dq  for  LaO  by  mass-effusion  than  by  mass-spectrometry  corresponds 
to  a  possible  presence  of  the  dioxide  which  was  not  corrected. 

Ackermann  et_  al.   (1964)  had  studied  the  vaporization  of  Y2O3 
by  Knudsen  effusion.     The  corrections  which  correspond  to  assumptions 
1,  2,  and  3  in  table  10  had  negligible  effect  on  calculating  the 
dissociation  energy  for  YO. 

Semenov  (1965)  had  used  mass-spectrometry  to  measure  the  Sc0+ 
ion-intensity  which  corresponds  to  the  vaporization  of  SC2O3.     The  2nd 
Law  value  of  Dq  for  ScO  in  table  12  was  calculated,  using  eq.  27.  The 
3rd  Law  value  of  D°  was  calculated  from  the  K_  of  reaction  3,  using 
Semenov' s  values  or  the  pressures  of  ScO  and 0  which  he  had  calculated 
from  his  measurements  of  the  effusion  rate  for  SC2O3  and  which  he  had 
corrected  with  his  estimate  of  8  for  the  ScO/Sc  pressure  ratio. 

Alcock  and  Peleg  (1967)  had  measured  rates  of  effusion  for  Gd203 
and  Y2O3  by  torsion-effusion  and  L an gmuir- torsion  experiments.  They 
had  assumed  that  vaporization  occurs  as  reaction  3  and  had  plotted 
log(P/atm)  versus  1/T.     Their  "total  pressure"  was  interpreted  as 
the  sum  of  the  LO  and  0  partial  pressures  so  that  the  Kp  of  reaction  3 
was  evaluated  by  reading  their  small  graph. 
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Table  12.     Dissociation  Energies  Recalculated  from  Published  Weight-Loss  and  Mass  Spectrometric  uata 


Sample 


T/K 
Range 


Source 


DQ(LO)/kcal  mol" 
Znd  Law  3rd 


2400-2570 
2551-2567 
2555-2693. 


Mass.  Spect.   (10  dV) 
Effusion  Rate  (3  runs 
Knudsen  (5  runs) 


.} 


Semenov  (1965) 
Ames  et.  al.  (1967) 


172 


160.4 


159.4 
160.9 


YO 


2509-2720 
2250-2390 
2492-2697 
2220-2620 


Knudsen  (11  runs) 
Langmuir-torsion,  Knudsen 
Knudsen  (14  runs) 
Mass-Spect.  (11  eV) 


Ackermann  et.  al_.  (1964) 
Alcock  et.  al_.  (1967) 
Ames  e_t.  al.  (1967) 
Ackermann  et.  al.  (1973) 


185 
173 
181 
171 


167.2 

162 

170.5 


LaO 

La  +  La203 
La203  +  Ta 


La2Oj 

La203  +  Re 


1650- 
1980- 
2234- 
2015- 
1490- 
2150- 
2258- 
1778- 


1900 
2220 
2441 
2280 
2190 
2560 
2427 
2308 


Mass  Spect.  (40  eV) 
Knudsen 

Knudsen  (22  runs) 
Mass  Spect.   (25  eV) 
Mass  Spect. 
Transpiration 
Knudsen  (5  runs) 
Mass-Spect.  (11  eV) 


Chupka  et.  al.  (1956) 
Kulvarskaia  et.  al.  (1960) 

Goldstein  et.  al.  (1961) 

Bushkovich  et.  al.  (1967) 
Benezech  et.  al.  (1969) 
Ackermann  et.  al.  (1971)^ 


259 
200 
187 
181 
194 


191.1 

197 

192.7 


199 
192.1 


2060-2400 


Kulvarskaia  et .  al.  (1960) 


214 


180 


NdO 

Nd203  +  Ta 


2080-2670 
2255-2434 
2095-2520 
2210-2560 


Knudsen 

Knudsen  (30  runs) 
Mass  Spect.  (25  eV) 
Transpiration 


} 


Kulvarskaia  et.  al.  (1960) 
Goldstein  et.  al.  (1961) 
Benezech  et.  al.  (1969) 


205 


166 
168 


169 
166 


178 


Sm203 
Sm203 


W 


2120-2350 
2333-2499 
2410-2670 


Knudsen 

Knudsen  (9  runs) 
Transpiration 


Kulvarskaia  et.  al_.  (1960) 
Ames.  e_t.  al.  (1967) 
Benezech  et.  al.  (1969) 


179 

138.9 

130 


143 

138.1 

133 


Eu203 
Eu203 


Ta 
W 


2050-2300 
1984-2188 
2192 


Knudsen 

Knudsen  (9  runs) 
Knudsen 


Kulvarskaia  et.  al.  (1960) 
Ames  et.  al.  (1967) 
Petzel  et.  al ,  (1971) 


139 
94 


125 

131.6 

131a 


GdO 
Gd203 


+  Ta 

+  W 


2080-2380 
2160-2390 
2408-2546 
2350-2590 
2380-2670 


Knudsen 

La ngmui  r- torsion 
Knudsen  (9  runs) 
Knudsen  (18  runs) 
Transpiration 


Kulvarskaia  e_t.  al.  (1960) 
Alcock  et._al.  (1967) 
Ames  et.  al.  (1967) 
Messier  (1967) 
Benezech  et.  al.  (1969) 


222 

114 

183 

171.4 

169 


177 

157 

169.2 

165.8 

172 


DyO 
Dy  0 
Dy203 


+  Ta 
+  W 


2260-2460 
2432-2637 


Knudsen 

Knudsen  (16  runs) 


Kulvarskaia 
Ames  et.  al 


et.  al.  (1960) 
(1967) 


16 
158 


155 
148.1 


HoO 

Ho203 

Ho203 


+  Ta 
+  W 


2230-2490 
2487-2711 


Knudsen 

Knudsen  (13  runs) 


Kulvarskaia 
Ames  et.  al 


et.  aj..  (1960) 
,  (1967) 


131 
48 


156 
153.1 


Er20-  +  Ta 
Er203  +  W 


2270-2490 
2492-2687 


Knudsen 

Knudsen  (16  runs) 


Kulvarskaia 
Ames  et.  al 


et.  al.  (i960) 
(1967) 


154 
170 


157 
151.4 


2450-2641 


Knudsen  (11  runs) 


Ames  et.  al.  (1967) 


163 


134.8 


Yb  03 
Yb203 


+  Ta 
+  W 


2060-2400 
2371-2626 


Knudsen 

Knudsen  (12  runs) 


Kulvarskaia 
Ames  et.  al 


et._al.  (1960) 
(1967) 


144 
127 


123 

93 


LU2U3 

Lu.2-3 


+  Ta 
+  U 


2120-2400 
2615-2700 


Knudsen 

Knudsen  (9  runs) 


Kulvarakaia 
Ames  et.  al 


et.  al.  (1960) 
■  (1967)  


84 
128 


183 
157.2 


As  reported  by  source 
Reference  2. 
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Ames  et  al_ .   (1967)  had  measured  the  rates  of  effusion  for  several 
sesquioxldes .     Except  for  YbO  the  values  of  D~  in  table  12  are  the 
results  of  assumption  3  described  for  table  10.     In  the  case  of  YbO, 
assumption  5  described  for  table  10  was  used. 

Bushkovich  et  al.   (1967)  had  tested  their  mass-spectrometer  by 
subliming  La203.     They  had  measured  the  LaO+  ion-intensity  as  a 
function  of  temperature  and  deduced  a  AH^  to  be  4  35  kcal  mol~l, 
presumedly  for  reaction  3. 

Messier  (1967)  had  measured  rates  of  effusion  for  Gd20.j  in  a 
tungsten  cell  with  0.47  cm  radius  for  which  the  Clausing  factor  was 
estimated  to  be  0.53.     This  cell  was  fitted  with  either  of  two  lids 
with  orifices  which  had  radii  of  0.450  and  0.855  mm  and  Clausing 
factors  of  0.356  and  0.507,  respectively.     Using  the  average  rate 
of  effusion  through  each  orifice,  simultaneous  solution  of  eqs.  10 
and  11  gives  the  vaporization  coefficient  cx  as  0.074  and  the  vapor 
saturation  _f  for  each  cell  as  0.96  and  0.82,  respectively.     The  D° 
for  GgO  in  table  12  was  calculated  by  using  assumption  3  described 
above  for  table  10. 

Benezech  and  Foex  (1969)  had  used  the  entrainment  method  and  a 
rotating  solar  furnace  to  study  the  vaporization  of  sesquioxides  which 
served  as  their  own  containers.     Thev  had  assumed  that  vaporization 
occurs  as  reaction  3  and  had  calculated  pressures  of  the  monoxides  of 
La,  Nd,  Sm,  Dy,  Tb ,  Gd,  Er  t  ard  Ho,  which  they  gave  on  a  small  graph. 
At  this  time  no  attempt  was  made  to  read  the  graph;  however,  the  authors 
did  report  2nd  and  3rd  Law  values  of  the  AHq  of  reaction  3,  involving 
La,  Nd,  Sm,  and  Gd,  although  they  gave  no  sources  for  the  values  of 
the  Gibbs  free  energy  functions  they  had  used. 

Petzel  and  Greis  (1971)  had  used  the  Knudsen  effusion  method 
to  study  the  influence  of  tungsten  cn  the  vaporization  behavior  of 
Eu203  at  2192  K.     They  had  placed  be  th  W  and  Eu203  in  an  iriduim-lined 
cell  to  control  solid-solid  contact,  had  analysed  the  effusate  for 
a  Eu/W  molar  ratio  of  3.7  to  3.8,  and  had  calculated  partial  pressures 
of  various  tungsten  oxides  and  the  cissociation  energy  of  EuO  which 
is  given  in  table  12. 

Recently  Ackermann  and  Pauh  (refs.  2,  3,  and  4)  have  used  both 
Knudsen  effusion  and  mass-spectrometry  to  study  the  vaporization  of 
some  sesquioxides  both  alone  and  in  the  presence  of  their  metals.  They 
had  also  studied  several  isomolecular  oxygen  exchange  reactions.  The 
values  of  D°  for  LaO  in  table  12  were  calculated  from  their  rate  of 
effusion  and  ion-intensity  data  which  corresponds  to  the  La203  vaporized 
in  a  Re  cell.     Only  corrections  for  the  metal  atom  and  diatomic 
oxygen  were  made.    A  review  of  ref.   3  was  not  completed  in  time  for 
inclusion  here.     In  ref.  4  Ackermann  and  Rauh  had  used  only  mass- 
spectrometry  to  study  the  vaporization  of  Y2O3.     Their  YO  ion- 
intensity  data  was  used  to  calculate  from  eq.  27  the  2nd  Law  value 
of  D°  for  YO  in  table  12. 
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Chapter  10 


Bibliography  on  Spectroscopy  of 

Fluorides  and  Oxides 
Belonging  to  Lantbanide  Series 

This  bibliography  lists  papers  cn  the  spectroscopy  of  lanthanide 
fluorides  and  oxides  in  gas,  liquid,  or  solid  phases.     It  also  includes 
some  papers  dealint  with  the  chemiluininescence  or  quenching  of  these 
compounds  in  excited  states. 

The  list  is  divided  into  fiftenc  sections: 

I .  Cerium 

a.  Cerium  fluorides 

b.  Cerium  oxides 

c.  Cerium  oxide  fluorides 

II.     Dyspros  ium 

a.  Dysprosium  fluorides 

b.  Dysprosium  oxides 

c.  Dysprosium  oxide  fluorides 

III.  Erbium 

a.  Erbium  fluorides 

b.  Erbium  oxides 

c.  Erbium  oxide  fluorides 

IV.  Europium 

a.  Europium  fluorides 

b .  Europium  oxides 

c.  Europium  oxide  fluoride; 

V.  Gadolinium 

a.  Gadolinium  fluorides 

b.  Gadolinium  oxides 

c.  Gadolinium  oxide  fluorides 

VI.  Holmium 

a.  Holmium  fluorides 

b.  Holmium  oxides 

c.  Holmium  oxide  fluorides 
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VI I •  Lanthanum 


a*   Lanthanum  fluorides 

b.  Lanthanum  oxides 

c.  Lanthanum  oxide  fluorides 


Till,  Lutetlum 

a*  Lutetiun  fluorides 
b,    Lutetlum  oxides 


IX.  Neodymlum 

a*   Neodymlum  fluorides 

b.  Neodyaiua  oxides 

c,  Neodymlum  oxide  fluorides 


X.  Promethium 

a,  Promethium  fluorides 

b.  Promethium  oxides 


XI,  Praseodymium 

a.  Praseodymium  fluorides 

b.  Praseodymium  oxides 

c.  Praseodymium  oxide  fluorides 


XII,  Samarium 

a.  Samarium  fluorides 

b.  Sanarium  oxides 

c.  Samarium  oxide  fluorides 


XIII,  Terbium 

a.  Terbium  fluorides 

b.  Terbium  oxides 

c.  Terbium  oxide  fluorides 


XIV.  Thulium 

a.  Thulium  fluorides 

b.  Thulium  oxides 
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XT.  Ytterbium 


a*  Ytterbium  fluorides 
b.   Ytterbium  oxides 


The  list  1b  based  on  the  files  of  the  Chemical  Kinetics  Information 
Center  and  an  examination  of  Chemical  Abstracts*  (1967-1974) • 

Suggestions,  additions  and  corrections  are  welcome. 


The  Chemical  Kinetics  Information  Center,  Institute 
•    for  Materials  Research,  National  Bureau  of  Standards, 
Washington,  D.C.    20234,  is  ^art  of  the  National  Standard 
Reference  Data  System.    The  Center  prepares  bibliographies 
for  general  distribution  and  also  in  answer  to  specific 
requests. 

The  Center  does  not,  normally,  distribute  copies  of 
published  papers  or  reports,  although  it  will  assist  in  the 
location  of  obscure  or  hard-to-obtain  material.  Papers 
should  be  obtained  locally  or  from  the  author  or  publisher. 
Gorernment  R&D  reports  that  are  available  to  the  public 
may  be  purchased  from  the  National  Technical  Information 
Service,  Springfield,  Virginia  22151. 


The  Center  indexes  and  stores  copies  of  published  papers 
and  reports  of  research  on  chemical  kinetics,  photochemistry 
and  collision  dynamics.  Contributions  to  its  collection  are 
welcome. 
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I,  Cerium 


(a).    Cerium  fluorides 

Batsanov,  S.  S#,  Derbeneva,  S.  S.,  and  Batsanova,  L.  R.t  Electronic  Spectra 
of  Rare  Earth  Metal  Fluorides,  Qxyfluorides  and  Oxides,  Zh.  Prikl. 
Spektrosk.  10.,  332  (1969)  (Russ);  Chea.  Abstr.  70:101332k  (1969) 

Baunan,  R.  P.,  and  Porto,  S.  P.  S.,  Lattice  Vibrations  and  Structure  of 
Rare-Earth  Fluorides,  Phys.  Rev.  161,  8*2  (1967);  Chea,  Abstr. 

67: 112549a  (1967) 

Claudel,  J.,  Hadni,  A.,  and  Striraer,  P.,  Low-Frequency  Vibration  of 
Lanthanide  Fluoride  Lattices,  C.  R.  Acad.  Sci«,  Ser.  B  274.  943 
(1972)  (Fr);  Chea.  Abstr,  77:40868k  (1972) 

Hadni,  A.,  Striaer,  P.,  and  Vermillard,  F,,  Far  Infrared  Absorption 
Observation  of  the  Low  Frequency  Lines  of  Praseodymium  Fluoride 
and  Other  Rare-Earth  Fluorides,  C.  R.  Acad,  Sci.,  Ser.  3  gg, 
1113  (1969)  (Fr);  Chem.  Abstr.  72:49216a  (1970) 

Ovsyannikova,  I.  A.,  and  Nasonova,  L.  I,,  Ljjj  X-Ray  Absorption  Spectra 

of  Praseodymium  and  Cerium  and  Effective  Coordination  Charges, 

Zh.  Strukt.  Khinu  XL,  548  (1970)  (Russ);  Chem.  Abstr.  73:82137u  (1970) 

Wertheim,  G,  K#,  Cohen,  R.  L. ,  and  Rosencwalg,  A.,  Multiplet  Splitting 
and  Two-Electron  Excitation  in  the  Trivalent  Rare  Earths,  Electron 
Spectrosc,  Proc.  Int.  Conf.  1971,  p.  813  (1972),  Shirley,  D,  A,,  Ed., 
Chem.  Abstr.  77:158530x  (1972) 


(b).    Cerium  oxides 
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